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ABSTRACT 


We  discuss  integrable  Hamiltonian  systems  -of  the  form; 

n0  n j 

V 


H(x,y)  = j Tj  + V V(x  -x  ) + V W(x  ) , 
i=l  l<i<)<n  ' 


(V(x),  W(x))  of  the  form: 

2 
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(A)  (x  ”,  xZ),  or  (B)  ((-  coth  | )2,  a-  eX)  , 

( 

-^and  also  systems  whose  Hamiltonians  are  integrals  of  these  systems^The 

derivation  of  «w  results  is  made  possible  by  the  fact  that  the  equations  of 
motion  can  be  interpreted  as  deformation  equations  for  matrix  functions  whose 
spectrum  remains  fixed  as  the  system  evolves  in  time,  leading  both  to  inte- 
grals of  the  motion,  and  a description  of  the  solution.  In  particular  we  show 
that  for  case  (A),  a purely  imaginary,  all  solution s/^^he  flow  are  periodic, 
addition,  for  all  these  systems,  symmetric  polynomials  in  the  x^(t)  turn  out 
to  be  rational  in  exponentials,  or  polynomials  in  at  most  n exponentials.  We 
then  study  the  scattering  maps  associated  with  these  systems  for  a > 0, 
obtaining  striking  properties  for  these  maps  and  interrelations  between  the 
systems. 
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SOME  TINITE  DIMENSIONAL  INTEGRABLE  SYSTEMS 
AND  THEIR  SCATTERING  BEHAVIOR* 

Mark  Adler 


1.  Introduction. 

We  conside"  Hamiltonian  systems  of  n particles  on  a line  interacting 
with  each  other  where  the  Hamiltonian  is  of  the  form: 


(1.1) 


1 n 2 n 
H (x,  y)  = - V y.  + V V(x.-x  )+  Y W( x , ) . 

" 2 i = 1 1 l<W<n  1 J itl  1 


The  examples  of  such  pairs  of  potentials  (V(x),  W(x))  to  be  considered  are: 

.?  /■ 2 
(A) 


(x  , - — x ) 


a 2. 

T 


(B) 


((jcoth  |)2,  a eX  ) . 


Calogero,  Marchioro^  and  Sutherland^  have  studied  some  of  these 

potentials  in  the  context  of  quantum  mechanics,  and  their  work  suggested 

looking  at  the  classical  systems.  For  the  case  a = o 
[31 

J.  Moser  has  shown  that  both  of  the  above  examples  are  integrable 
systems,  i.  e.  , possess  n integrals  whose  associated  Hamiltonian  flows 
commute,  and  in  addition  the  integrals  are  rational  in  (x.,  y.),  (eXi,  y.) 
respectively.  The  method  he  used  was  based  on  the  iso'spectral  technique 
of  Lax^  first  applied  by  H.  Flaschka^  ^ to  the  Toda  lattice.  This  consists 
in  the  construction  of  a matrix  function  of  (x,  y)  whose  spectrum  remains 
fixed  in  t if  x = x(t),  y = y(t)  are  solutions  of  the  above  Hamiltonian 
system.  We  then  take  the  eigenvalues  of  the  matrix  to  be  the  desired  integral 
and  study  systems  whose  Hamiltonians  are  functions  of  these  integrals. 
Sponsored  by 

1)  the  United  States  Army  under  Contract  No.  DAAG29 -7  5-C -0024 

2)  th^>  National  Science  Foundation  under  Grant  No.  MCS75-17  385 

3)  the  Office  of  Naval  Research  under  Contract  No.  N00014-76-C-0  3-01. 


We  extend  this  method  to  some  new  systems.  Moreover,  we  con- 
struct a second  matrix  function  of  (x,  y,  ti  whose  spectrum  is  invariant 
under  the  Hamiltonian  flow,  which  allows  us  to  describe  the  solutions  more 
or  less  explicitly. 

In  this  way  we  show  for  solutions  x = x(t)  of  the  above  Hamiltonian 

system,  in  case  (A),  a + 0,  that  the  symmetric  homogeneous  polynomials 

in  x.(t)  of  degree  v are  polynomials  of  degree  v in  e*0*  [6]  , which 

implies  that  all  solutions  of  the  system,  for  a purely  imaginary,  are  periodic. 

Furthermore,  for  a class  of  systems  whose  Hamiltonians  are  integrals  of  the 

above  system,  the  symmetric  polynomials  of  x,(t)  are  polynomials  of  at 

most  n exponentials  eV  . The  same  result  is  true  for  system  (B),  for 

a = 0,  if  we  replace  x.  by  e*i,  as  first  proven  by  Olshanetzky  and 
$ $ 

Perelomov  , and  moreover  for  systems  which  we  will  construct  such  that 

their  Hamiltonians  are  integrals  of  (B),  a = 0 . While  for  (B),  a t 0,  we 

find  that  the  symmetric  polynomials  in  e i are  rational  in  n exponentials 
X t 

e 1 , . . . , e n . 

It  is  then  easy  to  discuss  the  scattering  behavior  of  the  above  systems 
quite  explicitly  in  case  the  particles  ultimately  disperse,  as  in  the  case 
a-  > 0,  and  to  construct  scattering  maps.  For  instance  in  case  (A),  a > 0 
the  solutions  behave  asymptotically  like; 


(1.2)  x (t)  = -±-  {q*e±at  + «-1pS:F£>t  + 0(e-2,Q't| 
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for  t -*  ± or,  and  for  all  k,  q*  < q±  if  x (0)  < x (0)  . We  include 

K Krl  K K + l 

-1  -1 

the  artificial  factors  2 z,  a 2 ?,  for  a later  purpose,  to  make  (x,  y)  -*  (p,  q) 

a canonical  transformation.  On  the  other  hand  for  case  (B),  <y  > 0,  the  solutions 
behave  asymptotically  like; 

(1.  3)  xv(t)  = ± \fct  + pi  + 0(t"S,  for  t -*  ± oo  , 

where  \.  < X_  < . . . < \ <0,  if  x.(0)  < x . for  all  i . 

12  n ’ l i+l 

The  scattering  maps  which  we  construct  are  canonical,  given  by  poly- 
nomial relations,  and  lead  to  surprising  algebraic  transformations  between  the 
above  systems  (1.  1).  Moreover,  they  are  found  to  agree  wdth  their  own  inverse, 
i.  e.  , they  are  involutions.  For  instance  in  formula  (1.  3)  we  find 

(1.  4)  Pk  + = 2 log(«  . 

It  is  surprising  that  the  scattering  map  for  the  system  (A),  with  a > 0 , 
which  relates  data  at  t = - oo  with  data  at  t = + or  is  precisely  equal  to  the 
scattering  map  for  the  same  system  (A),  with  a - 0,  but  which  relates  data 
at  t = 0 with  data  at  t = » . Another  interesting  fact  is  that  the  scattering 
map  for  the  system  (A),  with  a > 0,  which  relates  data  at  t = 0 with  data 
at  t =30,  transforms  the  Hamiltonian  of  that  system  into  an  integral  of  the 
system  (B)  with  a - 0,  after  a trivial  change  of  coordinates. 

In  Section  IV  we  show  that  Hamiltonian  systems  near  by  to  system  (1.  1), 
case  (A)  with  a purely  imaginary,  have  at  least  n geometrically  distinct 
orbits  on  each  energy  surface.  While  such  results  are  easily  derived  in  case 
the  periodic  solutions  of  the  unperturbed  system  are  isolated,  they  are  quite 
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delicate  for  manifolds  of  periodic  orbits,  which  we  encounter.  This  result 
follows  directly  from  the  perturbation  theory  of  A.  Weinstein  [8]. 

I wish  to  thank  Professor  J.  Moser  for  his  many  suggestions  in  both 
formulation  and  technical  specifics,  and  last  but  not  least  for  his  enthusiasm 
and  encouragement. 


* Part  of  this  research  was  done  while  at  the  Courant  Institute  of 
Mathematical  Sciences  [14],  N.  Y.  U.  , with  the  Office  of  Naval  Research, 
Contract  No.  N00014-76-C-03-01. 

**  I am  indebted  to  F.  Calogero  for  communicating  to  me  the  then  un- 
published  results  of  Olshanetzky  and  Perelomov  at  the  above  mentioned  con- 
ference. Their  work,  which  was  done  independently  has  meanwhile  been 
published  in  [11].  Part  of  their  results  overlap  with  some  results  presented  in 
Sections  2-4  of  this  paper.  The  proofs  presented  are  those  of  the  unpublished 
preprint  [12],  which  was  handed  out  at  the  above  mentioned  conference  of 
January  1976,  in  conjunction  with  a brief  research  announcement. 

***  This  result  was  motivated  by  a recent  note  by  Sawada  and  Kotera.  See  [9]. 
****  Some  of  this  work  was  announced  and  presented  at  the  conference  on 
Theory  and  Application  of  Solitons,  held  January  1976  in  Tuscon,  Arizona,  and 
will  appear  in  the  proceedings  of  this  conference  [13]. 


-4- 


2.  Integrability. 


To  construct  the  integrals  of  H^(x,  y)  defined  by  (1.1),  Case  (A), 
we  use  the  matrices  of  J.  Moser's  system,  Y = Y(x,  y),  B = B(x),  where 
(x,  y)  « ST  = {(x,  y)  c R2n | x.  < x for  i = 1,  2,  . . . , n-1 },  S2  defined  once 

and  for  all. 

V = [L(*>y)1jk  = 6ikyi  * 1(1-V  ‘W’1 

V‘V  <XJ-Xk>'2)  > 


and  the  diagonal  matrices  X = diag(x  , x , . . . , x ),  D(y)  = diag(y  , y , . . . , y ) 
For  future  use,  we  observe  the  crucial  fact  that  L(x,  y)  has  n distinct  real 
eigenvalues.  (See  [3]). 

We  observe  the  commutator  relation  [X,  Y]  = C,  C = i(l  -6  ),  noting 

J K j K 

that  iC  is  the  identity  operator  on  the  subspace  in  Rn  specified  by 


i=l 


n 


Define  T (G  9 - G 9 ) to  be  the  Hamiltonian  vector  field 

G ft  yi  xi  xi  n 


acting  on  functions  of  (x,  y),  i.e.  if  x =G  , Y.  = -G  , i = l,  ...  n , 

1 yj  1 xi 

y)  ^ * 

then  — mr2 = OF  . As  an  operator  on  matrices  D_  acts  componentwise. 

at  (j  G 

Let  6 be  the  operator  acting  on  matrix  functions  of  (x,  y),  and  t,  defined 
by 


6(Z)  = (Dr  Z - f B,  Z])  + 9tZ  . 
a 

Note  that  6 is  a derivation,  i.  e.  , it  satisfies  the  sum  and  product  rule  of 
differentiation,  respecting  order,  for  it  is  the  sum  of  three  derivations.  We 
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observe  that  in  the  case  where  t doesn't  occur  explicitly  in  Z,  then 

6Z  = D Z - [B,  Z]  . We  also  note  for  future  reference,  that  I shall  stand 
H 

a 

for  the  n X n identity  matrix. 

9H  9H 

Theorem  1.  If  x,  y obey  x.  = - — - , y,  = - * i = 1,  . . . , n,  then  the 

1 dy , l ox , 

matrices  M*  = (Y  ± »X)  e* a satisfy  the  isospectral  differential  equation 

(2.  1)  6M±  = 0 , 

and  consequently  so  does  the  time  independent  matrix 
E = j M+  M'  = j (Y  + a X)  (Y  - aX)  . 

Proof.  We  establish  the  following  partial  differential  equation. 

(2.  2)  a)  6X  = Y,  b)  6Y  = aZ  X . 

Note  if  we  neglect  the  commutators,  these  equations  become  scalar  harmonic 
oscillator  equations.  Since  x = y,  we  have  X = D(y),  and  (2.  2)a  is  just 

D(y)  = Y + [B,  X],  i.e.  , Y = D(y)  + [X,  B],  and  since 

[*-  BV  = lx)  - V V = “l  - V 'VV'1  ■ 

the  definition  of  Y yields  the  result.  To  prove  (2.  2)b,  we  note 

H = l-[tr  Y2  - a2  tr  X2]  = H f P , 
a L U 

and  so  D Y = D Y + D Y,  but  in  [3]  it  is  computed  that  D Y = [B,  Y] 
a 0 0 

(a  straightforward  computation),  and  since 

2 

DpY  = Dp(D(y))  = a X,  we  have 
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Dpj  Y - Y f Dp  Y = [B,  Y]  + o"  X,  and  (2.  2)b  is  proven, 
o 0 

Now  (2.  2)  immediately  implies 

6(Y  ± aX)  = ± o(Y  ± oX)  , 

hence 

6M±  = (6(Y  ± aX))  e*01  f (Y  ± oX)(6(e T at  I)) 

= ± cr(Y  ± aX)  e Tat  + (Y  ± oX)(T  ae?at)  = 0 , 
and  thus  SM*  = 0 . 

That  E satisfies  6E  = 0 is  an  immediate  consequence  of  6 being  a deriva- 
tion. 

Remark  1. 

We  not  that  the  operator  equation  6D=  0,  or  equivalently 
U (6D)  U = 0 = ^-(U  ^DU),  where  U = BU,  is  equivalent  to 

U(t)  1 D(t)  U(t)  = U \o)D(0)  U (0),  provided  U(t),  U(t)'1  exists.  We  say  L undergoes 
an  isospectral  deformation  in  time.  Thus  the  constants  of  the  derivation  6,  i.  e.  , the 
algebra  which  forms  the  kernel  of  the  operator  6,  are  nothing  but  the  iso- 
spectral matrices  which  evolve  through  their  similarity  class  via  the  infinitesimal 
generator  B . Their  eigenvalues  and  all  functions  of  them  remain  constant  in 
time. 

Rema-k  2. 

❖ 

In  our  case  B f B =0,  so  apriori  U(t)  exists  for  all  finite  time,  and 
UU  = I if  it  holds  at  some  time  t . If  in  addition  U(t)  converges  for 
t -*»,  we  may  set  the  value  of  U at  t = *>  to  be  I,  the  n xn  identity 
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matrix.  Thus  U (t)  = I - f B(x(s))  U(s)ds,  and  we  now  redefine  U(x,  y)  as 

t X 

the  solution  of  the  integral  equation  U(x,  y)  = I - J B(x(s))  U(x(s),  y(s))ds  , 

‘ 0 

as  a function  of  x = x(0),  y = y( 0)  on  ft,  rather  than  along  the  orbit.  In 

other  words  we  construct  U(0)  via  the  above  integral  equation,  thus  directly 

exhibiting  its  dependence  on  ft  . Then  U(x,  y)  is  a solution  of  the  partial 

differential  equation  D U = BU  . 

H 


a 

Thus  for  every  point  (x,  y)  of  ft  we  have  U(x,  y),  B(x),  and  that 

U \x,  y)  L(x,  y)  U(x,y)  is  constant  along  orbits  and  equals  lim  L(x(t),y(t))  . 

t -*  so 

If,  as  in  Case  (A),  a = 0,  lim  L(x(t),  y(t) ) is  a diagonal  matrix,  then  the 

t -*  x 

columns  of  U(x,  y)  are  the  eigenvectors  of  L(x,  y)  . 

If  in  some  way  we  can  find  a normalization  of  the  eigenvectors,  and 

hence  of  U(x,  y),  we  may  dispense  with  its  definition  using  the  partial 

differential  equation  in  (x,  y)  space.  For  instance,  if  B satisfies  an 

equation  of  the  form  CB  = 0,  with  C a constant  matrix  with  no  column 

identically  zero,  then  automatically  GU  = C is  the  desired  normalization. 

For  = D (CU)  = CBU  = 0,  but  lim  U(x(t),  y(t))  = I by  definition, 

d t H 

t—  30 

and  Cl  = C,  hence  CU(x,  y)  = C for  all  (x,  y),  i.  e.  , U(x,  y)  is  a member 

of  the  Lie  subgroup  of  the  unitary  group  specified  by  CU  = C . This  symmetry 

coupled  with  the  fact  that  U diagonalizes  L,  algebraically  defines  the 

matrix  function  U if  the  eigenvalues  are  distinct.  The  above  happy  situation 

of  simple  eigenvalues  occurs  for  all  cases  we  will  consider,  and  C will  be 

of  the  form  (C  ) = i,  L shall  have  simple  spectrum,  and  so  we  shall  always 
K J 

be  able  to  define  our  U by  algebraic  processes. 
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Remark  3. 


We  note  that  E does  not  contain  time  explicitly  and  so  I (E)  = tr  EJ  , 
j = 1,  n are  n rational  integrals  of  the  motion. 

Letting  z = (x,  y)  , we  define 

vzg)=  £;  ||L2L-  , 

i.e.  , the  Poisson  bracket  of  f,  g,  computed  in  z = (x,  y)  coordinates, 

where  ( , ) is  the  ordinary  scalar  dot  product  in  R^n,  the  gradient 

operator  in  R“n  equipped  with  coordinates  z = (x,  y),  x,  y e R , and  J the 

2n  X 2n  constant  matrix  [;:•]  , with  I the  n X n identity  matrix.  Note 

D f = {f.g}  by  the  previous  definition  of  D . 
g g 

Theorem  2. 

The  I . { E)  are  in  involution,  i.  e.  , (I  , I } = 0 for  j,  k = 1,  2,  . . . , n . 

3 3 k 

The  I are  algebraically  independent,  i.e.,  the  system  H in  Case  (A)  is 
classically  integrable. 

To  prove  the  theorem,  we  use  the  asymptotic  description  (a  > 0)  of 

the  orbits.  Clearly  it  suffices  to  prove  that  4> , , = {I.,  I.  } vanishes  in  some 

3k  j k 

neighborhood,  since  6 (x,  y)  are  rational  functions  of  their  arguments. 

J K 

Also  it  suffices  to  consider  the  case  of  a > 0 , as  a enters  rationally  into 

* . The  algebraic  independence  of  the  I (x,y)  is  easily  seen  upon  letting 

J K J 

x.^j  - x.  be  large  for  all  i . For  a > 0 we  use  the  asymptotic  behavior  of 
the  solution,  letting  a = 1 for  simplicity. 
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Lemma  2.  1. 


For  (q,  < q < . . . < q ),  (p  , p , . . . , p ) varying  in  some  appropriate 
12  n 1 2 n 

neighborhood  N of  R2n,  there  exists  solutions  of  our  system  for  t > 0 , 
satisfying: 


T Xj(t,  p,  q)  = -p-  (q^  e1  + p^  e_t)  + 0(e~2‘  St)"| 

|Z-  3'  ) I t -t  -2.5.  f for  > 

y.it,  p,  q)  = — (q.  e - p.  e ) . °<e  ) J 

k 1 . / o 1 1 


where  the  above  expression  may  be  differentiated  with  respect  to  ( p,  q>  and 
still  remains  valid.  Moreover,  the  map  t:  (p,  q)  -*  (x(0),  y( 0))  is  canonical. 
We  sketch  the  proof.  For  u.(t,  p,  q)  = x.(t,  p,  q)  - — - (q.  e*  + p e ) , 

3 J ^7  3 3 

we  have  the  following  incegral  equation 


(2.4) 


u.  = J K(s-t)  V.fu,  s]ds  , 
1 *-  1 


where  K(t)  = sinht  , 

9H0  1 t -t 

V.fu,  s]  = --r— ^ (u(s,  p,  q)  + (q  e tpe  ))  , 

3 9Xj  si  2 

u(t,  p,  q);  R + X N - W — Rn, 

where  Rf  = { s | s > 0},  and  N is  the  neighborhood  in  R2n  to  be  determined. 

For  e > 0,  we  apply  the  contraction  principle  to  the  set  S = 9(e)  of  vector 

C2  functions  u = u (t,  p,  o)  of  W,  with  the  restriction  that  sup  ( I u | + 

W 

|i^|)e2‘St=  ||u||  < e • We  then  pick  M sufficiently  large  depending  on 
9t  — £ 

e,  fix  a q°  having  the  property  mini  a”  - qcj|  > M£,  q°  < q^,  and  define 

1 ^ 

N={(q  p)l  | p | + | q - q °|  < M }.  Then  one  verifies  in  a standard  manner 

t 

that  KfBjC  f B],  and  IlKfu^  - K[uj||  < 1/2  lluj-uJI  for  Uj,  u2  e g . 
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Hence  via  the  contraction  principle,  u = Kfu]  possesses  a unique  solution 


9 U 

u . If  we  define  the  vector  v = (u,  — . — l 


9q  ’ 9p”  = S^P(ivl  + l|f,)e2'5t 

then  a similar  argument  with  this  new  norm  will  yield  the  existence  of  — — 

3q 

9u 

— , their  derivatives  with  respect  to  time,  and  the  asymptotic  estimates 


on  these  quantities  stated  in  (2.  3). 


Thus  a solution  to  (2.  3)  has  been  constructed,  and  only  the  last 
statement  of  the  Lemma  2.  1 remains  to  be  proven.  We  now  have  the  map 
t , where 


, x(t,q,  P). 
y(t,q,  PK’ 


t >°, 


(q,  P)  « N . 


By  the  statement  of  the  first  part  of  the  lemma,  for  t sufficiently  large, 
t is  C'  invertible,  on  perhaps  a smaller  neighborhood,  and 


(2.  5) 


p),  y(t,g,  p)> 
9(q,  p) 


Now  the  map  (x(0),  y(0>)  -*  (x(t),  y(t))  is  canonical  and  hence  C* 
invertible,  thus  tq  = <(>  * is  C ' invertible  in  N . If  we  denote  the 
Jacobian  matrix  of  these  maps  by  j,  we  have 

a)  j(TQ)  = j (4>  S • j(Tt)  , 

T -2.  5t 

b)  j (Tt)I  i(Tt)  = -I  + 0(e  ),  as  a consequence  of  (2.  5), 

T t t t 

c)  j (4>  )J  j(4>  ) = J,  since  4>  is  a canonical  map, 

d)  llj^SlIlCe111  . 

T - 5 1 1 1 

As  a consequence  of  a)-d)  we  conclude  J(TQ)  J j(TQ)  = -J  + 0(e"*  ) for 
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T 

all  t . Hence  j(TQ)  J -Htq)  = ~J>  anc^  thus  T = T0  *s  canonical,  and 
the  lemma  is  proven.  We  observe  a map  preserves  the  form  of  the  previously 
defined  Poisson  bracket,  if  and  only  if  its  Jacobian  at  every  point  belongs 
to  the  symplectic  group  Sp(2n,  R)  of  2n  X 2n  matrices  satisfying  A JA  = J , 
which  is  easily  verified  using  the  definition  of  {,  } . More  generally  we 
consider  a transformation  whose  Jacobian  A satisfies  A JA  = pj,  with  a 
constant  p * 0,  and  call  it  a canonical  transformation  with  multiplier  p . 

Proof  of  Theorem  2. 

So  far  we  have  proven  the  map  (x(0),  y(0))  - (q,  p)  is  canonical,  but 
with  multiplier  -1,  i.  e.  {,}=-{,},  r = (q,  p)  . We  define  the  map 
4-:  (q,  P)  - (£,  q)  by 

(2.  6)  = log  q^  r][  = -qi  p.  , 

in  the  domain  0<  q^  < . . . < q^  . This  map  is  clearly  canonical  with  multi- 
plier -1,  and  hence  iJj  °t  \ t(N)  -*  4j(N)  is  canonical  with  multiplier  1 . 
Now  we  recall  that  (a  = 1,  for  simplicity), 

E = M +M"  = j(Y  + X)  (Y  - X)  = ^-(Y2  - X2)  + l-  [X,  Y]  , 

hence 


2Eu*6k/(yk-xk  + £'(xk-  v'2> + 

m 


. V v-,'  1 

(1  - (i:  :'x  ) + 1 + 1 'xk - XJ  (x,-xm> 

K I m 


-1 
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I 


where  the  prime  indicates  that  undefined  terms  are  to  be  left  out.  Using 
(2.3),  (2.6),  we  conclude 


(2-  7) 


. i ^ i 

lim  E(x(t),  y( t ) ) = EM  = (6k/  + -(1  - 6^)  coth 5 + 7O  - 6UJ)} 


= z(r,,e)  + j c = z (T!,e), 


2'  k r 


s,  = i<‘  - v 


where  Z{£,,r\),  z (it  h)>  are  defined  by  this  relation,  and  we  define 


z (n>  i)  - z(o>  t)  - 1C  ' Now  ^ 


(2.  8) 


F<i.  e - 1 tr  - j 1 \ + j l <=°th2 


i.  e.  , F(r|,  i)  - H^(£,  ri)  of  (1.  1)  case  (B) , a = 0,  then  the  flow  defined  by 


’i  dF 


'i  9F ' 


ds  9ri . i’  ds  9$,  * 
'1 


i = 1,  • • • , n , 


can  also  be  expressed  by  an  isospectral  deformation  of  the  form 


Z = D Z = [K,  Z],  see  [3].  Here  K has  off  diagonal  elements 


(2.9) 


Kk,U)  • Kj  slnh-'t^-e,)]2  = i qk  q(tqk  - q,f2 


and  the  sum  of  the  rows,  and  columns  is  zero,  as  was  the  case  for  B . 
Hence  CK  = -iK  = KC  implies  fC,  K]  = 0 = D^C,  and  therefore  Z+(t),£)  = 
Z(tj,  £)+jC,  (and  similarly  Z~(t],  £)  = Z(tj,  £)  - j C),  satisfies  the  same 
isospectral  differential  equation  as  Z . In  addition,  under  the  flow  induced 
by  F,  we  find 

i{  = Q[  s + + 0(  s*1)  , 


hi  = ai  + °(  s ) , 


(2.  10) 


for  s 


+ oc 


t 


with  a.  < o , and  in  fact  the  transformation  (£,,  n)  - (P,  °)  is  canonical 

by  a similar  argument  as  in  Lemma  2.1,  see  [3]  for  details.  However,  we 

will  sketch  a proof  that  the  o.'s  are  in  involution.  By  (2.7),  (2.10),  the 

distinct  oj  s constitute  the  simple  spectrum  of  the  matrix  Z+(r),  £)  (and 

also  Z ( r| , 4 ) ) » which  undergoes  an  isospectral  deformation.  Hence 

£>  = (£,r|),  is  a constant  of  the  motion  by  the  Jacobi  identity, 

which  we  evaluate  at  s=®,  But  since  n.  — a , for  s -*  oo,  in  the  C' 

i i ’ 

sense,  0 = {r|.,  <*j  and  thus  0 = {<*.,  a ^ } . Now  since  the 

matrices  E(x,  y)  and  lim  E(x(t),  y(t))  = Z+(q,  £ ) and  lim  Z+(n( s),  £(s))  are 

t-*^  S-^ori 

unitarily  equivalent  via  remark  1,  while  all  the  maps  considered  are  canonical, 

we  conclude  that  the  eigenvalues  of  E(x,  y),  which  are  the  o, a , are 

1 n 

in  involution,  and  hence  so  are  the  I (E)  = V a being  functions  of  the 

j k=l  k 

eigenvalues.  Thus  the  proof  of  Theorem  2 is  completed. 


* 
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3.  Scattering  Theory. 


To  discuss  the  scattering  theory  of  the  system  of  Section  2,  we  need 

a lemma,  which  will  be  a consequence  of  the  following*. 

dH  dH 

Theorem  3.  The  solutions  of  the  system  x.  = —2-  y = - - — - 

l 9y.  i 9x.  ’ 

l l 

i = 1,  . . . , n,  obey 
(3.  1) 


n 


y (2s-v)art 


' x.  (t)  = ) c e' 

H , l u „ s 
1=1  s=0 

( 

v a natural  number,  c a rational  function  of  initial  data.  Hence  the 
’ s 

i Oft 

x,(t)  are  explicitly  given  as  algebraic  functions  of  e , and  initial 
data  (a  * 0)  . 


Proof.  Referring  back  to  the  definitions  of  Section  2,  one  computes 

,,  a at  -at. 

( 3.  2)  X = — — (M  e -Me  ) , 


since  M*  = (Y  ± aX)  eTat 


th 


Raising  both  sides  of  Equation  (32)  to  the  v power,  and  then 
taking  the  trace,  we  conclude 


(3.1) 


n v , , , 

V XV  = V c<v>  e(2s-v)a 

i — / ~ s ’ 

i=l  s=0 


( v)  V ( v) 

where  c = tr  C , and  C is  a sum  containing  all  possible  v-fold 
s s s 

products  of  the  noncommuting  matrices  Mf,  M taken  s,  and  v-s  times 
respectively.  We  recall  that  the  algebra  generated  by  Mf,  M undergoes 
an  isospectral  deformation  during  the  flow,  by  remark  1 of  Section  2.  Thus 

( v] 

C undergoes  an  isospectral  deformation  during  the  flow,  and  so 
s 
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(t)  = c*l\o)  is  an  integral  which  is  a rational  function  of  the  initial 

s s 

data.  We  amplify  Lemma  2.  1,  and  now  prove  a statement  about  the  scattering 

of  all  solutions.  Since  the  particles  can't  collide  we  may  order  them  so 

that  x,  < x.  , for  i = 1,2 n-1.  and  all  real  t. 

li+l  ’ ’ ’ 

Lemma  3.  1.  For  a > 0,  every  solution  satisfies 


[3.  3) 


...  1 at  -1  -at  - 

xk(t)  = 3—  (qk  e - a pk  e ) + 0(e 

•si  2 


2at 


) , 


as  t — +oo,  < dj  + j,  j = 1,  2,  , n-1,  under  the  hypothesis  of  Theorem  3, 


-1 


with  a similar  conclusion  for  t-»  -»  . 

Proof.  We  note  in  Theorem  3 that  c^  = tr  (Yn  + aXl]  , where  the 
subscript  indicates  evaluation  at  t = 0 . Since  Y ± aX  = L(x,  y ± ax)  , 
we  conclude  Y ± aX  has  simple  spectrum,  for  we  observed  at  the  beginning 
of  Section  2 that  matrices  of  the  form  L(x,  y)  have  simple  spectrum.  Multi- 


-avt 


plying  (2.  1)  by  e , 
(3.4) 


we  find 


-at.  v (v)  v -2 at  -4at. 

> (x.  e ) = c ' + c , e + 0(e  ) , 

j v v-1 


with  c 


= (*/T)'  V V q V,  where  qt  < q^  < . . . < qn 


are  the  distinct 


i 


,-l 


time-independent  eigenvalues  of  M . Since  x,  x2,  ...,x  are 

2 n v 

distinct  and  ordered,  they  are  uniquely  determined  by  ^ x , and  hence 

<v,  ' = ’ 1 

by  the  expression  for  c^  , (3.4),  and  by  the  distinctness  of  the  x.'s 

and  q 's  respectively,  we  have 
i 


Xie 


’at  = — (q.  + a_1p.  e"ZQ,t)  + 0(e"4at)  . 

■s/T  1 1 
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This  defines  the  p.  . Upon  multiplication  of  both  sides  of  the  latter  ex- 
pression by  eat,  the  lemma  follows. 

Remark  1. 

We  thus  have  a map  <j>  : 12  -»•  12,  12  defined  in  the  beginning  of 

Section  2,  given  by  4>  + (x,  y)  = (q,  p)  for  a real.  We  wish  to  show  it  is 
a bisection  of  12  . Given  q.  < q , Pj,  we  may  form 

ao 

x.(t)  = (q.e0,1  + a'  p e‘at  + f sinh(s-t)  V.[x(s)]ds,  i = 1,  2,  . . . n , 

n/T  1 t 1 

and  prove  the  existence  of  a solution  x ( t) , for  t > tQ,  with  the  given 

(q,  p),  by  the  methods  of  Lemma  2.  1.  From  the  conservation  of  the  energy, 

n 2 

H (x,  y),  and  the  behavior  in  time  of  given  by  (3.1),  we  conclude 

n°  I t | i=l 

Yj  ( I Xi  I + I y.  I ) £ C e a , and  therefore  x(t)  exists  and  is  unique  for  all 
i=l  1 

time,  and  in  particular  at  t = 0 . Whereas  Lemma  3.  1 enables  us  to  define 
<i>+  on  12,  the  above  discussion  shows  that  is  well-defined  on  12,  and 

hence  4>+  Is  a bijection  of  12,  and  canonical,  by  the  arguments  of  Lemma  2.  1. 

We  now  exhibit  an  important  feature  of  4>^,  namely  that  it  is  an 
algebraic  mapping,  i.  e.  , it  is  implicitly  given  by  polynomial  relations, 
and  it  is  an  involution,  i.  e.  , 4>+  = identity.  Then  we  shall  discuss  a 

surprising  relation  of  4>^  with  the  scattering  map  in  the  case  a = 0,  and 
rational  invariants  of  this  map. 

The  algebraic  nature  of  the  map  is  given  by  matrix  equivalences. 

Referring  to  remark  2 in  Section  2,  since  |b|  = 0(t’  ) along  an  orbit,  B 
being  the  infinitesimal  generator  of  the  isospectral  flow,  we  may  define  the 
function  IJ(x,  y),  and  by  (3.  “i)  and  the  same  remark,  in  a purely  algebraic 
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manner;  and  thus  conclude  lim  M*  (t)  exists  and  equals  U \x,  yjM*  (x,  y)U(x,  y). 


Since 


t-*  or 


M*  = (Y  ± aX)  eTat  , 


x ryf  —1  X at 

[M  (z(t),t)]kj  = 6kj(yj<t)  ± ox.(t))e  + i(l-6jk)(xk(t)  - x.(t))  e 


and  upon  using  y (t ) = x.(t),  (3.  3),  we  find  for  t -+<*>, 

(y.  + o x.)(t)  e'at  \IZ  a q , ( y . - o-x^tje"1  - - \I  2 pA  , 


e at(xi<t)  - x.(t))'1  -*  0, 


- Xj)(t)  -*•  \1  2 (c?i  - ^j)1  J » 


and  so  we  have  (see  remark  2 of  Section  2), 


(+)  U *(x,  y)  M+(x,  y)  U(x,  y)  = lim  M+  (z(t),  t)  = 2 qQ 

(3.5)^  _ t-oe  _ _ 

(-)  U (x,  y)  M (x,  y)  U(x,  y)  = lim  M-(z(t),t)  = - V 2 P } , 

t-»00 


where  Q = Diag  (q  , . . . , q^),  P = L(q,  p)  . 

It  is  important  to  notice  that  both  matrices  M±(x,  y)  are  transformed 

---IT 

in  (3.  5)  by  the  same  similarity  transformation.  Now  let  L (x,  y)  = 2 2 a L (x,  y) 
-I  -1 

= 2 p a L(-x,  y),  and  let  ~ denote  spectral  equivalence,  remembering 
X = diag(Xj,  . . . , x^),  we  obtain: 

Theorem  4. 

The  scattering  map  4>f  is  given  by  the  following  symmetric  relations 
( 3.  6a)  X ~ L (q,  p f «q) 

a 

( 3.  6b)  L^(x  y f a x)  ~ 0 , 

a 

where  both  equivalences  are  effected  by  the  same  similarity  transformation 
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T -I 

U (x,  y)  . From  the  symmetry  of  (3.  6)  we  see  at  once  = <j>^,  i.  e.  , 

4>  o 4>  = id  . 


Proof.  Line  (3.  6b)  is  the  transpose  of  equation  (3.  6)^,  while  ( 3.  6a)  is 
the  transpose  of  the  difference  of  equations  (3.  6)*,  thus  finishing  the 
proof. 


We  now  study  the  scattering  map  4>  from  t -*•  -00  to  t -*•  4oc, 

<j>  = » (<t>  )-1,  where  4>  denotes  the  map  from  the  initial  data  (x(0),  y(0)) 

into  the  scattering  data  (q  , p ) for  t = -so  , defined  by 
x = _J_  (q‘  e‘at  + a'lp*  e^S  + 0(e2a,;),  t - -»  . To  compute  we 

1 n/T  1 1 

may  restrict  ourselves  to  the  case  a - 1 . We  recall  that  Lj(x,  y)  = 

2 ~ ? L (x,  y),  and  we  conclude  from  Theorem  4,  that  <j>+>  4>  _ are  given 


implicitly  by : 


(3.7) 


X = U + L1(q+,  p+  4 n + )U~*  , 
L^(x,  y 4 x)  = U+  Q+  U"\ 


(b) 


X = U_  L^q",  P‘ 
L^x,  -y  4 x)  = U 


q'u;1  J 


We  may  derive  (b)  from  (a)  by  using  the  time  reversibility  of  the 
differential  eauations  under  (t,  x,  y)  -*  (-t,  x,  -y),  hence  <}>  = <j>^  ° p , 

where  p(x,  y)  = (x,  -y)  . 

Therefore  we  conclude  from  (3.7), 


U Q"  U-1  = L^(x,  -y  + x)  = L^(x,  -(y  + x))  + V2  X 

-r  T -1  +T  — T + 4 + . T 

= -L.(x,y  + x)  + v/ 2 X = (U  + ) f-0  4 \!  2 (q  , P 4 a ) } U + 

= (uy1  {sir  Ly\ pf)>  y , 
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and  after  a similar  calculation  involving  Q+  we  conclude; 

Theorem  s.  For  case  (A),  a > 0,  the  scattering  map 

= <t>  + ® <t>  is  given  implicitly  by 

0’  ~ LT(q+,  P+),  L(q ",  p")  ~ Qf  , 

where  both  of  the  equivalences  are  brought  about  by  the  same  unitary  map, 

V = UI  U.  ’ 

The  case  a - 0 has  been  treated  by  Moser,  but  for  completeness 

we  include  short  proofs  in  the  spirit  of  the  above  arguments,  and  then 

relate  the  results  of  the  two  cases.  It  was  shown  in  [3],  that  if  x,  y 
9H0  9Ho 

°bey  *i  = ’ *i  = - . i = l,...,n,  and  if  x^O)  < x1+1(0),  i = 1,  . . . 

then  as  t oo 

( 3.  8)  Xj  = t f p.  + Oft"1),  y.  - q.  + Oft'2),  qt  < q^,  i = 1,  2,  . . . , n-1 

In  fact,  if  we  define  the  map  by  4>+(x,y)  = (q,  p),  4>+ : u - S2  where 
U is  the  domain  used  in  the  definition  of  <t>  , then  4,  is  seen  to  be  a 
canonical  bijeetion  by  arguments  similar  to  ones  given  here.  We  are  now 
in  a position  to  state; 

Theorem  f>.  = 4>,  where  <$>  = 4>+  • $ 1 . 

Moreover  4>  = <t>  f 0 = 4>  + ° p 0 4>  f so  ^ is  ‘conjugate1  to 

the  linear  reflection  p,  and  hence  is  an  involution. 

The  equation  4,  = expresses  that  the  scattering  map  4>  relating 

data  at  t = - 00  to  data  at  t = +<*  for  the  system  H , is  precisely  the 

a 
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scattering  map  relating  data  at  t = 0 to  data  at  t = » for  the  system 
H . Even  the  unitary  maps  in  the  two  cases,  which  effect  the  matrix 
transformations  are  the  same,  as  they  both  are  contained  in  the  Lie  subgroup 
discussed  in  remark  2 of  Section  2. 


Proof.  We  now  observe  that 

,+ 


-1 

a . . . + 


M = M (a),  - M ) = X coshart 


sinhat 


Y , 


are  isospectral  matrices  under  the  flow  Hff,  (see  remark  1 of  Section  2), 
and  for  « - 0 we  obtain  L(x,y),  X - t L(x,  y)  as  isospectral  matrices  under 
the  flow  induced  by  HQ  . Using  (3.8),  and  the  arguments  of  the  Lax  formalism 
used  previously,  we  conclude: 


(3.  9) 


lim  L(x,  y)  ~ L(x(0),  y(0)) 

t-*  on 


9 


} LT(q,  P)  = L(-q,  p)  = lim  (X  - tL(x,  y)) 
V t— V 


9 


where  ~ denotes  unitary  equivalence. 

Hence  (x(0),  y(0))  = (q,  p)  is  explicitly  given  by 
X(0)  ~ LT(q,  P),  L(x(0)  , y(0))  ~ 0 , 

where  both  equivalences  are  effected  by  the  same  similarity  transformation. 
This  completes  the  proof  of  the  theorem,  by  comparing  this  with  Theorem  5. 

The  map  linearizes  the  flow  for  the  case  a = 0,  and  is  a alge- 

braic, yet  we  have  only  an  implicit  description  of  it.  However,  we  can  find 
some  rational  functions  invariant  under  the  action  of  4j+,  similarly  for  <t>  + • 
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Theorem  7. 


(i)  Let  V be  the  algebra  of  real  valued  rational  functions  of  (x,  y)  e S2  , 
generated  by 

trff(X)  • f(L  (x,  y fax))],  (tr  f(X) ) • tr  f(L  (x,  y fax)))  , 

a a 

where  f is  a polynomial.  Then  every  member  of  V is  an  invariant  under 
the  action  of  <f>  . 

(ii)  Similarly  let  V be  the  algebra  of  real-valued  rational  functions 
generated  by 

tr[f(X)  • f(L(x,  y))],  (tr  f(X))  • (tr  f(L(x,  y)))  , 

f a polynomial.  Then  V is  a list  of  rational  invariants  of  4*  . 

Proof.  Let  be  the  algebra  of  matrix  valued  functions  M(x,  y)  of  (x,  y)  , 

T 

generated  by  X = diag(x  . . . , x ) and  L (x,  y + ax)  . Note  that  any  homo- 

In  a 

T 

morphism  of  is  completely  specified  by  its  behavior  on  X,  L (x,  y 4-  ax) 

a 

through  extension.  Now  the  map  4>+(x,  y)  = (q,  p)  induces  a mapping 

T 

(X,  L^(x,  y + ax))  - (L(q,  p 4-  a q),  0)  , 
explicitly  given  by  (3.7a).  By  the  above  remark,  this  extends  naturally  to  a 

A 

homomorphism  4>  : 17[  -*■  7!\ , which  by  (3.  7a)  is  expressed  by 
(3.10)  $+  (M(x,  y))  = (x,  y)  M(x,y)  ty(x,  y)  . 

Since  4>+  (x,  y)  = (q,  p),  we  can  consider  functions  of  (x,  y)  as  functions  of 
(q,  p),  which  we  shall  do  when  we  think  of  ^ as  the  range  of  the  map  <£+  , 
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i.e.,  we  shall  regard  $ (M(x,  y) ) as  a function  of  (q  , p)  . We  note  that 

A A 

trace  is  an  invariant  of  4>  . Hence  defining  f(x,  y)  for  any  polynomial 

f by  the  following  line: 

f(x,  y)  = tr[f(X)  • f(L^(x,  y + ax))]  = tr($  + ff(X)  • f(L^(x,  y + ax))]) 

T 

= trff($+(X))  • f($+(La(x,  y + ax)))]  = tr[f(L^(q,  p + aq))  ■ f(Q) ] 

= tr[f(Q)  • f(L  (q,  p + a q))]  = f (q,  P)  , 

a 

we  find  f (x,  y)  = f (q,  p)  if  (x,  y)  = (q,  p),  i.e.,  f(,  ) is  an  invariant. 

The  other  part  of  (i)  is  immediate,  while  (ii)  is  proven  precisely  as  (i)  . 

We  remark  that  an  alternate  derivation  of  Theorem  7 can  be  given, 
whereby  the  invariance  of  the  above  quantities  is  seen  to  follow  from  the 
'invariance'  of  the  Poisson  bracket  under  canonical  transformations,  (see  [14]). 
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4.  Periodicity  of  Solutions  for  a Purely  Imaginary  and  a Perturbation  Result. 


By  Theorem  3,  the  homogeneous  symmetric  polynomials  of  degree  v in 

£ at 

the  x.  are  polynomials  of  degree  v in  e , for  the  system  (A),  a 1 0 . 

On  the  other  hand,  the  x.  are  strictly  ordered,  x^  < x2  < . . . < xn»  due  to 
the  singularity  in  the  potential,  hence  the  x.  are  uniquely  given  as  algebraic 
functions  of  e±ot  (and  the  initial  data).  We  thus  conclude,  verifying  a 
conjecture  of  Calogero's  [10]: 

Corollary  3.  2. 

9H  9H 

If  (x,  y)  obey  x . = ~ , y . = - , i = 1,  . . . , n,  a pure  imaginary, 

then  the  solutions  are  all  periodic  with  (not  necessarily  primitive)  period 
2niO' 

Remark  1. 

Since  the  motion  is  periodic,  one  expects  2n-l  integrals  of  the  motion 
to  exist,  which  in  this  case  would  be  rational  in  (x,  y),  and  indeed  they  do. 
Take  for  the  first  n integrals  I^[E]  = tr  E , and  for  the  latter  n-1  integrals, 
Real[I  (M+)  1^  (M  ”)],  J = 2,  n . One  sees  they  are  algebraically  in- 
dependent by  considering  the  x far  apart,  in  which  case  the  matrices  are 
nearly  diagonal. 

Remark  2. 

A linearizing  transformation,  which  is  rational  in  one  direction,  alge- 

( ***)  , v ^n 

braic  in  the  other  direction,  is  given  by  r:  (x,  y)  — (tr(L  + «X)  j = 

(Zj,  . . . , z ) = z,  (x,  y)  « 52,  52  defined  at  the  beginning  of  Section  2,  and 
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maps  the  flow  into  z . = jaz,,  (see  [9]).  The  mapping  t may  not  be  in- 
vertible, but  for  x. - x,  large,  L + aX  is  close  to  diag(x^  + a y^,  . . . , x^  + oy^ ) 
and  thus  by  the  implicit  function  theorem,  knowledge  of  z enables 
us  to  recover  (x,  y)  from  z and  z . Hence  z , ...,z  are  algebraicly  in- 
dependent, complex  valued  rational  functions  of  x,  y . 

Theorem  8. 

2 

For  every  sufficiently  small  C perturbation  of  H , there  exists  at 

a 

least  n geometrically  distinct  periodic  solutions  on  every  energy  surface 

E > E = min  H (x,  y)  . 
a a 7 

(*,  y) 

While  such  results  are  easily  derived  in  case  the  periodic  solutions 
of  the  unperturbed  system  are  isolated,  they  are  quite  delicate  for  mani- 
folds of  periodic  orbits  which  we  have  encountered.  We  shall  use  Theorem  1.  i 
of  f8],  a beautiful  theorem  of  A.  Weinstein.  Instead  of  stating 
Theorem  1.4  in  its  full  generality,  which  requires  much  terminology  [8], 
we  shall  just  list  an  immediate  consequence  of  it,  which  will  suffice  for 
our  purpose,  and  then  we  shall  verify  that  the  necessary  hypothesis  are 
satisfied  in  our  situation. 

Cor,  of  Thm.  1.  4,  [8] 

( e) 

Given;  The  Hamiltonian  system  H = H + e P(x,  y,  e),  (x,  y)  in  a 

neighborhood  of  the  manifold  given  by  the  relation  H^  = E,  and  VH°  * 0 

on  H^  = E,  P being  CZ  in  its  arguments,  while  for  the  value  E,  H^  = E 

is  a manifold  which  is  homotopic  to  the  sphere,  free  of  equilibria,  of  only 

(el 

periodic  solutions.  Then,  for  small  e,  the  system  with  E = H contains 
at  least  n geometrically  distinct  periodic  orbits. 
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Proof  of  Theorem  8. 


We  shall  apply  the  corollary  to  the  case  H = H , a purely  imaginary, 

a 9 

E > E = min  H (x,  y)  . 

a a 

x>  Y 

Thus  we  just  have  to  check  for  E > E , the  manifold  H = E is 

a a 

diffeomorphic  (homotopic  would  suffice)  to  S2n_1,  with  no  fixed  points. 

Now 


1 rt  2 

H = i-  y y.  + V , 

a 2 — ' 1 a 


v = . y x2  + y — 1 

a 2 <->  i A x,  - 


i<i  (xi  ■ V 2 ’ 


and  we  see  that  V (x)  is  strictly  convex,  since 


(ti>  Vxx^  = “T  + tTl)  = -a  <Zhi  ) + 6 Z 

dt  '^i 


(hj  - -Hj)' 


(6i-e()4  >0  * 


for  purely  imaginary  a . Thus  V (x)  assumes  a strict  nondegenerate  minimum 

at  some  unique  x , where  V (x  ) = E . We  conclude  H (x,  y)  is  strictly 
a a a a a 

convex  in  (x,  y)  € ft,  and  now  consider  the  surface  H = E > E , which  by 

a a 

the  above  contains  no  fixed  points  of  the  flow.  The  smooth  surface  H = E 

a 

is  thus  the  boundary  of  the  convex  body  < E,  and  hence  is  diffeomorphic 

to  a sphere,  the  diffeomorphism  being  given  by  spherical  projection  through 

an  interior  point  of  the  body  to  some  large  fixed  sphere.  Clearly  VH  * 0 

a 

on  H = E,  for  E > E . 

a 


Remark  2. 

For  n = 2,  3,  one  can  show  that  for  energy  surfaces  near  the  equili- 
brium point,  one  has  distinct  orbits  of  primitive  periods  2iria"1  • {-  — -} 

1 ’ 2 ’ ' ’ ' ’ n ’ 
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and  Gallavotti  and  Marchioro  [7]  conjectured  this  to  be  true  for  all  n 

2 

Thus  H is  quite  different  from  the  system  H = — ^ y2  . — V v2 
° 2 ^ i Z Lt  i 

for  which  all  solutions  have  the  same  primitive  period. 
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5.  Associated  Flows. 


The  integrals  found  for  the  system  H , case  (A),  were  of  the  form 

a 

tr  EP,  (see  Theorem  2).  More  generally  = tr  f ( E) , for  any  polynomial 
f,  gives  rise  to  an  integral.  It  is  thus  natural  to  consider  H^  as  the 

Hamiltonian  of  a flow,  which  we  call  associated  with  H = tr  E . For 

a 

p 

these  more  general  flows  tr  E or  also  tr  f(E)  are,  of  course,  again  in- 
tegrals since  all  these  functions  are  in  involution,  and  therefore  E is 
isospectral  under  all  flows  . We  will  determine  another  isospectral 
matrix,  which,  however,  depends  explicitly  on  t . This  will  allow  us, 

as  in  the  case  of  the  system  H , to  describe  the  solutions  of  the  system 

a . 

explicitly.  The  result  is  stated  in  Theorem  10  below. 

We  will  determine  these  isospectral  matrices,  first  for  another  flow, 
and  later  translate  the  results  to  the  flow.  For  most  of  this  section  we 
will  study  this  other  system,  which  is  closely  related  to  the  Sutherland 
system,  (1.1),  case  (B),  a - 0,  and  carry  out  the  above  plan  for  it.  The 
system  to  be  considered  now  is  given  by  the  Hamiltonian: 


1 v z 

[5.  1)  G(q,  p)  = - Y (a  p ) 
1=1 


(q,  P)  « ft  , 


ft  defined  in  Section  2. 


Notice  that  F(q,  p)  is,  after  the  trivial  canonical  change  of  coordinates  (Z.  6), 

transformed  into  the  system  given  by  (2.  8),  i.  e.  H of  (1.  1)  case  (B),  « = 0 . 

a 

For  the  matrix  description  of  the  flow  we  introduce: 


(5.2) 


{R(q,  p))jk 


6 


jk(qj  Pj>  + z (1  - V 
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R = R k C where  C = j (1  - 6 ),  and  therefore 

jk  £ )k 


p) = 6jk  (t,j  pj* +i<i-6it) — 1 


ik  <y<y 


*\ 


(5.  3) 


R4",.(q,  p)  = 6jk(Pj  Pj)  + i 0 - 6llf) 


jk 


Jk  (q,  - q. ) 


’J  k j 

We  observe  that  K(q,  p)  = Z(|,  -q),  R*  - Z , with  Z,  Z defined  in  (2.7), 
and  (q  , p)  - (£,  n)  defined  in  (2.  6)  for  qi  > 0,  i = 1,  2,  . . . , n . Thus  by  the 
discussion  immediately  following  (2.10),  R1^  , for  the  case  q.  > 0, 


+ * 


i = 1,  2,  . . . , n,  has  simple  spectrum.  We  note  that  by  definition  (R  ) =R  , 
R+Q  = QR',  where  Q = diag(qj,  q , . . . , qn)  . Thus  it  is  clear  that 
Rf,  R"  have  the  same  spectrum  if  0 is  nonsingular,  or  equivalently  that 
tr(R+)V  = tr(R~)V,  for  all  v . Since  the  latter  equation  is  rational  in  (q,  p)  , 
it  holds  for  all  (q,  p)  « U . Similarly  Q R = R Q implies  Q f'(R  ) = f'(R  ) Q 
for  f a polynomial,  and  so 


(5.4)  qi[f'(R")]ij  = [f'(R+)]ij  qj  ' 

If  f is  real,  then  [Q  f'(R~)]*  = [f'(R  )]  Q*  = f(Rf)Q  = Qf‘(R  ),  i.  e.  , 

Q f'(R  ) as  well  as  f'(R+)Q  are  Hermitian,  and  so  in  particular  q.[f’(R  )].. 
as  well  as  q.[f‘(R+)]ii,  and  [f(R-)]. . = [f*(R+)]i.  are  real.  Thus  we  con- 
clude 

i Yj  tf'(R+)ljkRkj  = i{ff,(R+)R+ljj  - ff’(Rf)ljj  RjV 

= i {[g(Rf)]jj  - qjPj[f,(R+)]jj}  » 


with  g(s)  = sf  (s),is  purely  imaginary. 
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An  easy  calculation  shows  that  for  F(q,  p)  defined  in  (5.  1), 
G(q,  P)  = tr(y(R')Z  + cl)  = tr(j(Rf)Z  + cl),  c = - H . 

Theorem  9. 

it  , r,  • 8r  9F 

Under  the  flow  q.  = 7^-  , p^  = - - — , with  the  Hamiltonian 

1 i 

F = Ff  = tr  f( R *")  = tr  f(R_)  , 


which  is  associated  with  G(q,  p)  of  (5.1),  the  symmetric  polynomials  in 
q.  = dj(t)  are  rational  in  exponentials  e^j1,  j - 1,  . . . , n,  provided  q.  > 0 . 
This  is  a consequence  of  the  fact  that  the  four  matrices, 


(5.6) 


e-tf,(R+) 


Q, 


M" 


Q e 


— tf ' ( R ) 


f 


undergo  isospectral  deformations  as  (q,  p)  = (q(t),  p(t) ) evolve  under  the 
F flow. 

Proof.  These  matrices  satisfy  an  equation  of  the  form  L = [A,  L],  where 
A=A  = A , A=-A,  A defined  by  the  two  equivalent  relations 


Vq’ p) 


ff,(R'»jk  Rkj 


The  identity  of  A+,  A is  an  immediate  consequence  of  (5.4),  (5.  3),  while 

* f 

A = -A  follows  from  the  fact  that  f(R  )Q,  hence  Qf(R  ) is  Hermitian, 
and  from  (5.  5).  To  prove  the  matrices  of  (5.6)  are  isospectral,  it  is  suf- 


ficient to  prove 
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(5.  8) 


(a)  R f = [A-,  Rfl,  (b)  M"  = [A',  M"]  , 

± * 

because  taking  the  adjoint  of  (5.  8),  the  fact  A = A = -A  implies  a 
statement  analogous  to  (5.8)  for  R , M . We  first  prove  (a),  i.  e.  , 
DpR+  = [A  , R f ] . Using  (5.  3),  and  defining  f f ’ ( R )]  = t , we 
make  a few  preliminary  calculations: 


f Dr  qs a ®p  [trf(R“)]  a tr[f,(R")  ’ ap  R']  = qstss  "I 


(5.9)  < 


v qj(tjs'tsj) 

DpP5  = [trf(R  )]  = -p.t  + i T.  — t1 


F s . q 


s ss  H . .2  . 

J*s  (qg-q  ) J 


i(1-6jk} 

where  we  have  used  [3  R ] = 6Jk(P,6)s>  * ‘2  («Jj  \s  - \ 6js»  • 

S (qfV 


P,q, 


Hence  Dr  V * <»k  pk>'  = q k V qk  pk  = 1 A (tjk  ' V • 


= 1 n — -- 

j/k  (aj-qk)' 


8u.  (A',  »‘]kk=Zk  % ^k  - Rk,  V 


y 

LJ 


''Vkj\  f qj 


qk-qj I \qr%l  lqk-qj  qj-qk/jj 


qjtjk 


= 1 ^ ^ 2)  (tjk-tkJ^ 


f ~ 

and  so  we  have  shown  D^R^  = [A  , R ] , and  we  now  must  show 


DF  RJk  5 tA’’ RV  Mk- 

+ 1 qJ 

We  compute  DR,  = D 


F jk  F (qj-qk) 

Vu  q),qii))-qktkkM  ‘'vv'kk-v 

- 1 1 — — : - 3 I - 2.  » 


M)-qk> 


,qi  - V 


<VV 
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hence 


qi  Aik  _ 

lqj-qk 


- =[A,R  Ijj,,-  Dp  Rj  j,,  (defining  A by  this  equation) 


v (A"  R+.  - R.+  A",  ) - D_  R+. 
-J  js  sk  js  sk  F ]k 


V jk_ 

W 


s*j,  k 


(<k  - R)V + % - v 


1 q, 


(qj-V 


Vjs 

— — 

q 

s 

i — 

c r 

l 

<VV 

<w 

w 

'skqs 

<W 


lq]  qk 


(qj  - v2 


We  now  wish  to  show  A =0,  and  for  that  we  shall  use  the  addition  law, 

jk  ’ 


|qj  - V 


<qj  -qs’  ,qs-V  ,q)-qs)  <Vqk>  ’ 


inside  the  bracketed  term.  So  we  compute 


i)k  = v(<k-R)Vt,(Ajj-vtiqk<^^- 


1 l <(tjs  - lsk> 


s*j,  k 


'VV  <VV 


We  compute  the  bracketed  term  separately; 
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(,ls  ■ 'sk' 


Mrqs)  (qs-V 


t.  q t q to 
js  s sk  s js  s 

— + — — 


(skqs 

(q)-qs>  <Vqk>  <W  ' <W 


qi  ‘u 

IS  llj  -qs> 


- 1,  + 


Vsk 

<Vqk» 


t,  + >v? . 

Is  (qs-qk) 


Vsk 

(qj-qs) 


Ajs  " Ask  " l(t|s  Rsk  ■ RJs  'sk1  • 


Hence 


jk 


= -s  a;v  i . ",s  Rsk  - rj‘s  ,sk> 

s s s*j,k  1 


+ ljk(Rkk  “ RjJ)  + Rjk(tjj  ' tkk) 


= ‘fZ  Ai's  - X AsJ  + [f,<R')»  R'l 

s 1 s 


Jk 


= a;s-s  v = i£  aj's  ’ bvi5-7>  • 

s s s 


So  we  just  have  to  show  i ^ Ajs  = i.  e.  , that  A C = 

s 

for  all  j k . But  since  V A , = 0,  for  all  k,  we  have 

’ Lj  sk  ’ ’ 

- . * . ~ s 

0 = (C  A ) = A C,  as  was  to  be  shown,  concluding  the 

It  remains  to  prove  (5.  8b)  i.  e.  , M~  = [A_,  M’]  . 


0,  where  C „ 
jk 

C A =0,  and 
proof  of  (5.  8a). 


= i 
so 
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We  first  compute  (M  )’  = (Qe_t^  ^ * = 


(DFQ)e-tf‘'R')+0(e-tf'(R'>)-  = 


(DrQ)e'tf'(R  ] + Q{-f(R')e'tf'(R  > + [A-,  e'tf'(R)]}  , 
where  we  have  made  use  of  DpR"  = [A-,  R-],  which  implies  D g(R  ) = [A-,  g(R~)] 


On  the  other  hand, 

fA'.M']  = [A‘,Oe_tf,(R  ']  = QfA',  e’tf'(R  }]  + fA',Q]e 
and  thus  in  order  to  have  M = [A  ,M  ],  we  must  have 
(5.  10)  DfQ  = Qf(R')  + [A",  0]  . 

In  other  words,  we  must  have 


-tf'(R') 


(5.11) 


df  q)  = •)) 

[0,  *-|)k(‘-v  = (VV,)k(1-V=  q)  V'-V  ’ 


but  (5.11)  is  an  immediate  consequence  of  (5.9),  (5.7),  and  (5.3)  . We  have 
thus  proven  (5.  8),  and  the  statement  of  the  theorem  concerning  the  rational 
character  of  the  solutions  will  be  easily  shown  to  follow  from 


(5.12) 


Q(t)  = W(t)Q0  etf  (R0)  w'\t), 


where  W is  a matrix  satisfying  W^  = I,  W = A~W  , 


which  we  quickly  verify.  Indeed,  as  a consequence  of  (5.  8), 

M"  = Q e'tf'(R  ] = W Qq  W'1,  R'  = WRfl  W_1  , 


where  the  subscript  shall  now,  and  in  the  future  indicate  evaluation  at  t = 0 , 
and  no  subscript  indicates  evaluation  at  t . Hence 
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Q = WQq  W"1  etf'(R  ) = WQ0W'1(Wetf'(R0)  W_l)=WQ0  etf  (V  w'1  . 

**  12 

We  note  that  Olshanetsky  and  Perelomov  have  proven  (5.12)  for  f(s)  - J s 


in  slightly  different  coordinates. 


The  rational  character  of  the  solutions  now  follows,  for  by  previous 
remarks,  R*  has  simple  spectrum  < v2  < " • • < vn  for  0 < ql < q2  < ' * ‘ < Qn  ’ 
and  so  by  (5.  12), 


tr  Qv  = J q%  P (e 

i=l 


f'(v  )t 


. , f'(v  )t 

where  P is  a polynomial  of  degree  v in  the  n exponentials  e 
v 

Remark  1. 

In  the  language  of  remark  1 of  Section  2,  the  matrix  equations  of 
motion  under  the  F flow  may  be  expressed  as: 

(5.13)  (a)  6R*  = 0,  (b)  6Q  = Qf' (R~)  = f’(Rf)0,  (c)  6M±  = 0 , 


where  the  derivation  6 is  defined  by 


(5. 14) 


6(-)  - [A,  (* )]  + at(*)  ■ 


Remark  2.  If  we  introduce  x.,  y.  by  = e J , Pj  = -iy,  1 - • • • » n » 

a canonical  transformation,  then  iR  ^ becomes 


6 


ikyj 


x 

7<*  - 6jk'  co,(“ 


If  we  take  f(s)  = j sZ  in  Theorem  9,  then  all  solutions  are  bounded  as  a 
consequence  of 
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1_2  1 2 1 2 i j 

'I 'r  R = 2 1/1*7  5-  COt  1-2— : 1 

i = l i<j 

being  a constant  of  the  motion.  To  be  more  specific,  we  also  observe  that 

t R ” “ 1 

Rq  is  diagonalizable,  since  Q = WQ  e o W , with  W unitary, 
ix  t R ” 

Q = diag(e  1,  . . . , e n),  and  thus  the  matrix  e o is  bounded  for  all  real 

t R ” t R ” 

t . Write  e o in  diagonalized  form  and  observe  that  since  e 0 is 
bounded,  Rq  has  purely  imaginary,  not  necessarily  distinct  eigenvalues 
iVj,  . . . , iv^  . From  which  we  conclude,  using  (5.  12),  that 


Q ~ S e 


iD(v) 


l s’  e'Y  , 
j = l 


where  D(v)  = diag(v^,  . . . , v ),  and  the  S^'s  are  the  column  vectors  of 
S . Since  necessarily,  |det  s|  = 1 * 0,  we  must  have  * 0 for  all  j . 
Thus  all  solutions  of  the  flow 


VV  Vi=-Hx.’ 
I 1 


i - 1,  • . • , n , 


1 vn  2 1 vy  2 i 1 

h = 2,X/i  *7  E cot  <— > ■ 


are  quasi-periodic;  and  moreover,  in  special  cases  they  will  be  periodic, 
namely  if  the  v are  integer  multiples  of  a number  p # 0 . 

The  quasi-periodic  character  of  the  solutions  would  be  a consequence 
of  a well  known  theorem  of  Arnold  [6],  if  only  we  knew  that  the  gradients  of 
the  n Integrals  of  the  compact  system  H = j tr  R^  were  everywhere  linearly 
independent.  Since  we  don't  know  this,  it  is  conceivable  that  such  phenomena 


as  exceptional  points,  hyperbolic  tori,  etc.,  occur,  which  is  ruled  out  by  the 
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quasi-periodic  behavior  of  the  solutions.  Of  course,  the  rational  character 
of  the  solutions  is  a stronger  assertion. 


Corollary  1.  Under  the  same  hypothesis  as  Theorem  9,  with  f(s)  replaced 

f ± T 

by  f(-»s);  i.  e.  , H = trf(-»(R  ) ),  we  have  for  the  time  evolution  of 
(q,  p)  = (q(t),  p(t)),  with  the  Hamiltonian  Hf, 


L (q,  q + <*p)  = W{  — — {Q  e 

a slT  ° 


-atf'(-ffR')  -1  atf 1 ( -»R~)  . 

0'+a  e 0 L(qQ,  PQ)  }]W 


-1 


Proof . Since  QL(q,  p)  = R+,  Q = WQq  e'atf’  ('aR0)  W'1,  R+  = WR^  w'1 

imply  L(q,  p)  = W eatf  (_Q,R0)  L(q0P0)  W \ we  get  from  L^(q,  P + »q)  = 

— - — (Q  + a ^ L(q,  p) ) , by  substituting  the  above  expressions,  the  stated  result. 


Now  we  return  to  our  original  family  of  systems  H = tr(f(E(x,  y) ) , 
and  the  plan  of  the  beginning  of  this  section.  The  tool  is  simply  the  canonical 
transformation  used  to  prove  Theorem  Z,  which  takes  into  . 


Theorem  10. 

If  x = 9 Hf(x,  y),  y = -8  Hf(x,  y),  i = 1,  . . . , n , 
1 yi  1 Xi 

then  the  three  matrices 


E, 


e-atf'(E)(Y 


aX),  (Y  - arX)  e 


atf’(E) 


undergo  an  isospectral  deformation  with  the  same  unitary  generator  for  all 
the  deformations,  and  we  have  the  time  dependence  of  the  solutions  given 
by 


(S.  1M 


X(t) 


(Y 


<,X0)  ♦ (»X0-Y0)e“>f  <V) 


where 


indicates  unitary  equivalence. 
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, we  have  for  a > 0 , 


Proof.  By  (3.  5),  (3.  10),  and  M*  = (Y  ± aX)  e*  at 

(5.  16)  {Y  + aX,  -Y  + „X}  - K/T  «Q,  s[T  LT(q,  p)  } , 

4>  + discussed  briefly  in  the  proof  of  Theorem  7.  We  recall  that  acted 
on  the  ring  H of  matrix  valued  functions  of  the  variables  (x,  y)  = 4>+\q,  p)  , 
generated  by  Y + aX,  Y - aX,  by  extending  formula  (5.  16),  so  that  for  M(x,  y) « 7![  , 

A _ } 

we  have  4>+(M(x,  y))  = U (x,  y)M(x,  y)U(x,  y)  . As  a consequence  of  the  defi- 

A 

nition  of  4>+,  we  found  in  (3.6a) 

$ (X)  = L (q,  p f aq)  , 

and 

(5.17)  $+(E)  = -a(R‘)T,  where  E = j(Y  + aX)( Y - aX)  . 


We  also  note  that  the  mapping  4>f  is  canonical  with  multiplier  -1,  so  we 

replace  t by  -t  in  the  change  of  coordinates  (q,  p)  -►  (x,  y),  to  describe 

A -1 

evolution  in  time.  Now  we  apply  the  inverse  (<j>+)  of  the  above  transforma- 
tion to  derive  from  (5.16),  (5.17)  that 


(S/ho  e„.f'(-«R-,)Ti(Jt)-l(e«,f(-«(R-,T,0)  = 


ea(-t)f*(E)  _J_(Y  + aX)=  _L_  e‘atf'<E)(Y  + aX)  . 


But  by  Theorem  9,  0 eat^  ( ) iS  5^1(31-  to  Q , hence  also  to  — — (Y  + aX  ) . 

-atf'(E)  ^ 2a  ° ° 

This  shows  that  e ° (Y  + aX)  is  similar  to  (YQ  + aX^),  and  thus  is  an 

isospectral  matrix  for  the  flow  . Formula  (5.  17)  states  E is  similar  to 

- T - T 

-a(R  ) , which  by  Theorem  9 is  similar  to  -a(RQ)  , hence  E is  isospectral. 

Formula  (5.15)  then  follows  similarly  from  3.6(a)  and  Corollary  9.1. 
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At  first  Theorem  10  is  valid  only  for  a > 0,  but  as  all  the  relations 
involved  are  algebraic,  we  may  take  a arbitrary.  The  formula  for  X in 
Theorem  10  yields  the  explicit  solutions  of  the  motion,  and  as  in  Theorem  9, 
we  conclude 


tr(L  ± oX)1  = P*  <e*‘rtf',vs>)  , 

V 


for  H = H^,  (see  Theorem  9),  the  v^'s  distinct,  for  a > 0 . 

Corollary  10. 1.  If  H = Hf  = tr  f(E),  the  differential  equations  of  motion, 
written  in  matrix  form  are 

(5.18)  a)  6E  = 0,  b)  6X  = j [f(E),  X]  + \ [f  (E),  Y]  , 


with  [A,  D]+  = AD  + DA,  6(- ) = DR(- ) + 3t(*  ) - [B,  (• )],  where 


(5.19) 


a)  \ ---J  . 


n 


Jjk  (Xj-Xk)  Jk 


b)  Z B 


S = 1 


Sj 


Proof.  From  Theorem  10  we  have 
+ -ortf'(E) 


E,  M = e 


(Y  + arX),  M = (Y  - «X)  e 


atf'(E) 


are  isospectral  under  the  H flow,  with  one  and  the  same  unitary  generator. 
In  other  words, 

SM*  = 0,  6E  = 0 . 

for  some  B(t),  which  depends  on  initial  data,  which  for  the  moment  we  use 
in  the  above  definition  of  6 . On  the  other  hand,  6M*"  = 0 implies 
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bM4'  = e~aA  6(Y  + «X)  + e'"tf  (E)(-af'(E))(Y  + oX)  = 0 , 
from  which  we  conclude 

(5.  20)  6(Y  + aX)  = af'(E)  (Y  f aX) 

and  similarly  from  6M  = 0,  it  follows 
(5.21)  6(Y  - aX)  = -o(Y  - oX)  f'(E)  . 

From  (5.20),  (5.21),  we  compute 

6E  = 6[  j(Y  + aX)(Y  - aX)]=  j(6(Y  + aX))(  Y - oX)  + |<Y  + aX)(6(Y  - aX)) 

= f f*(E)  (Y  + «X)(Y  - aX)  - |(Y  + aX)(Y  - oX)  f'(E)  = o[f'(E),  E]  = 0 . 

Formula  (5.  19a)  follows  upon  subtraction  of  (5.  20)  from  (5.  21),  upon  writing 
out  the  full  expression  for  6 . 

We  need  only  show  that  we  can  impose  the  normalization  of  (5.  19b), 
which  as  we  saw  in  remark  2 of  Section  2,  which  we  shall  be  constantly 
referring  to,  is  equivalent  to  requiring  that  the  unitary  generator  of  the 
matrix  flow  is  contained  in  the  Lie  subgroup  of  the  unitary  group  specified  by 
CU  = C,  i.  e. , with  Lie  algebra  specified  by  CB  = 0 . From  the  proof  of 
Theorem  10,  it  follows  that  the  unitary  generator  of  the  Hf  flow  can  be 
written  as  the  product  of  the  unitary  generator  of  the  Ff  flow  of  Theorem  9, 
and  the  U’s  occurring  in  the  map  of  (5.  16),  as  displayed  in  ( 3.  10). 

By  remark  2,  the  latter  U's  (see  the  discussion  preceding  (3.  5))  are  con- 
tained in  the  above  Lie  group.  The  U of  the  Ff  flow  has  as  its  infinitesima 
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generator  the  A defined  in  (5.7),  which  clearly  satisfies  CA  = 0,  and  thus 
by  remark  2,  this  U also  is  contained  in  the  above  Lie  group.  Thus  the 
unitary  generator  of  the  flow  can  be  written  as  the  product  of  elements 
in  the  above  specified  Lie  group,  and  so  by  remark  2,  we  can  impose  the 
normalization  condition  (5.  19b).  We  also  observe  that  the  fact  B(t)  = 

B(x(t),  y(t)),  i.  e.  that  B is  really  a function  on  £?,  was  not  assumed 
a priori. 


Remark  3. 

For  H = = tr  f(Y),  a - 0,  (5.  18)  reduces  to 

(5.  22)  6Y  =0,  6X  = f'(Y)  , 


where  6(-  ) = D^( • ) + 9^- ) - [B,  (•  )]  • 

We  also  note  that  (5.  22),  (5.  19)  imply 


(5.23) 


6Y  = 0,  6(X  - tf'(Y))  = 0 , 


[f'(Y)] 


(5.24) 


Ik 


Jk  (Vxk)’ 

y b.  = o 
Y js 


j * k , 


. 


and  so  by  (5.23),  the  solutions  of  the  H flow  obey 

X(t)  ~ (XQ  + tf'(Y0)),  hence  V x|'(t)  = tr(XQ  + tf,(YQ))1',  w = 1, 


the  latter  equation  being  a generalization  of  a result  found  in  [9]. 
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6 . Another  Integrable  Sy stem 


In  this  section  we  discuss  another  integrable  system,  namely  (1.1) 
case  (B),  a * 0,  formally  analogous  to  case  (A),  but  physically  behaving 
like  a system  of  interacting  particles  under  the  influence  of  an  additional 
force  acting  from  the  right.  We  shall  make  this  more  precise  in  the  discus- 
sion of  the  scattering  theory  of  the  system.  Moreover,  we  shall  show  the 
solutions  of  the  system  to  be  rational  in  exponentials,  by  getting  a fairly 
explicit  formula  for  the  time  evolution  of  the  system. 

Theorem  11.  The  Hamiltonian  system 

, n , , , x -x  n x. 

(6. 1)  H (x,  y)  = - V y + j V sinh"  + a v e 1 

i=l  l<i<j<n  i=l 

9H  9H 

is  integrable.  If  x , = , y = - , i = 1,  . . . , n,  then  the  matrix 

yi  i 

1 ■f’  — XX  x x 

(6.  2)  W = — Z (X,  y)Z  (x,  y)  + ae  , e = diag(e  1,  . . . , e n) 

undergoes  an  isospectral  deformation 

D W = [K,  W],  where  Z*  = Z ± j C , 
a 

with  Z(x,  y),  C,  K(x)  as  defined  in  (2. 7),  (2.9).  Then  the  n algebraically 
independent  functions,  t.  = tr  W^,  j = 1,  . . . , n,  rational  in  (y.,  eXi)  , 
are  in  involution. 

Proof.  We  shall  actually  prove  this  result  for  the  system 

I ^ 21 

F (Q,  P)  = T V (Q.P.)  +T  J (7— 71)  + a q , (q,  p)  € i2  , 

a 2 1=1  11  l<i<J<n 
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where  q is  not  necessarily  positive,  and  then  perform  the  usual  canonical 
transformation, 

(2.6)  q.  = e 1,  q.P.  = y^,  in  the  component  of  £2  with  q.  ^ 0 
for  all  i . We  then  have 


(6.3)  Z±(x,  y)  = R±(q,  p),  eX  = 0,  K(x)  = A(q)  , 

with  R,  Q,  A of  Theorem  9,  f( s)  = j sZ  . We  compute  the  matrix  differ- 
ential equations  of  motion,  using  the  derivation  &(•  ) = (•  ) - [A,  (•  )]  , 


namely: 


a 


(6.4) 


(a)  6Q  = — (QR  + R 0)  > 

(b)  6R±  = -aQ  . 


Equation  (6.4a)  is  the  same  equation  you  would  get  for  a = 0,  namely 

equation  (5.  13b),  for  f(s)  = j s2,  since  FQ  = — tr(R  ) , OR  = R Q , 

and  D q = 3 F is  independent  of  a . Since  D = D + a ^ 3 , 

?a  1 Pi  a Fa  F0  i Pi 

while  the  off  diagonal  elements  of  R*  are  independent  of  p.,  equation 
(6.4b)  follows  from  Theorem  9 , in  particular  (5.13a),  in  the  case 
f(s)  = j s2  . Hence,  by  (6.4), 

6(  — R + R'  + oQ)  = j (6  R + )R"  + ~ R+  & R“  + obQ 

= -§QR‘  -f  RfQ  + |(QR"  * = 0 . 

and  so  6W  =0,  i.  e.  , D W = f A,  W]  . The  only  thing  left  to  prove,  re- 

° i 

turning  to  (x,  y)  coordinates,  is  that  the  (trW  )’s  are  in  involution,  as 
their  algebraic  independence  is  easily  seen  by  considering  the  case,  for 


-4  3 - 


r 


all  i,  of  y.  very  large,  ' xj  j very  lar9e»  and  > 0 very  large. 
As  in  the  proof  of  Theorem  2,  due  to  the  rational  character  of  the  integrals 
trW*  in  (y,  eXi),  it  suffices  to  prove  their  involutive  character  in  some 
open  neighborhood.  We  shall  pick  the  neighborhood  in  the  (x,  y)  space 
where  for  a > 0,  the  solutions  have  the  following  asymptotic  behavior; 

xi = Y + pi + \ < °’  Xi  < Xi+i  for  a11  1 * 

(6.  5) 

Yi  = Y + °(t"  t - +»  . 


As  in  Lemma  2.1,  such  a neighborhood  is  easily  found  by  converting  the 
differential  equations  of  motion  into  an  integral  equation.  Then  for 


(6.6) 


X - X 

ZJk  = V)  ti(1-Vooth<J2J1) 


we  have 


lim  Z(t) 
t — oo 


1 Yjk + to)k} 


„)k  = -i  » k >i  , 
= +1  if  k < j , 
= 0 if  k = j , 


X 

while  limit  e = 0 . Thus  the  definition,  (6.2),  of  W implies  lim  W(t) 
t — oo  t—  or 

exists  and  is  a function  of  the  \.'s  . We  then  apply  the  same  argument  as 

in  Theorem  2 to  verify  that  the  tr  W^‘ s are  in  involution.  Also,  as  in 

Theorem  2,  since  the  integrals  are  rational  in  a,  the  statement  of  involution 

for  all  a follows  from  the  statement  for  a > 0 . 


-44- 


We  now  explore  this  system  (6.  1)  in  a sequence  of  corollaries.  Since 


no  collisions  can  occur  we  order  the  particles  x.  < x.+^  for  i - 1,  2,  . . . , n- 
all  t . 

Corollary  11. 1 . AH  solutions  of  the  system  Hq,  for  o > 0,  behave 
asymptotically  in  the  following  manner; 

x.  = ± Xjt  + P*  + 0(t  S, 

(6.7)  t - ± \1  < \2  < . . . < \n  < 0 , 

yt  = * \ + 0(t'2)  , 


(6.  8) 


p*  + p~  = 2 log(a\.  ) 


In  words,  the  scattering  behavior  of  the  system  is  that  of  n de- 
coupled particles,  each  interacting  with  the  origin  under  the  potential 
eX  . The  effect  of  the  interaction  is  felt  through  the  strict  ordering  of  the 
terminal  velocities. 


Proof.  With  no  loss  of  generality,  assume  a = 1 . Then  dropping  the 

subscript  H = H , for  a = 1,  (given  by  (6. 1)),  we  have 
a 


x 


i 


9H 

a — = Y„ 
dy.  i 


x 


8H 

av 


i = 1,  . . . , n 


im  plies 

^ x — 2 

(6.9)  x = V 4>(x  -X  ) - e , 4>(x)  = -(coth  j)  (sinh  -) 

)*i  ' 

We  first  show  that  for  all  i,  lim  y,(t)  exists  and  is  nonpositive, 

t -►  x 

for  which  we  use  only  the  following  properties  of  $ : 


4>'(x)  < 0,  = -<t>(x),  x4'(x)  >0,  x f 0,  <(>(0)  = » 


From  (6.  9)  we  have 


*1  = - e 1 * 
J>1  ] 


and 


we  note  4>(Xj-x^),  -e  1 are  < 0,  while  the  energy  relation,  (6.1)  , 


gives  an  upper  bound  on  |x  .1  . Hence 


2 si 


2H  > | f x l dt  | = -J'  x dt  = V f <j)(x  .xp  dt  + J' 


e 1 dt 


J>1 


and  so  4>(x  -Xj)  > 0,  e 1 are  integrable  on  (-*>,  °°)  . From  this  it  follows 

Xj(<»)  = J x^(t)dt  + Xj(0)  exists.  Now  from  (6.9)  we  have  x ^ + 4>(Xj-x^)  = 

X x 

V 4>(x  -x  ) - e 2f  where  again  we  note  4>(x  -x  ),  -e  2 are  < 0 for 

J>2  2 1 

j > 2 . By  the  previous  step  and  the  upper  bound  on  fx^  |,  we  have  that 
the  left  side  of  the  above  equation  is  integrable  in  t on  (-«>,  <*>),  and  thus 
so  is  the  right  side.  Since  all  terms  on  the  right  hand  side  of  the  equation 
are  of  the  same  sign,  each  one  of  them  is  integrable,  i.  e.  4>(x^-x^)  , 

J * 2,  eX2  are  thus  integrable,  and  so  as  before,  x2(°°)  exists.  In- 
ductively we  climb  our  way  up  ic  i = n,  and  conclude  x . (°°)  exists  for 


all 


j,  and  eXj,  4>(x,-x.)  are  integrable  on  (-<», »)  for  all  i * j . Setting 

^ J nr 


\f  = X («-),  since  x.  < x and  f eXj^d t exists  for  all  j,  we  must 


J J' 


J J+l 


have  \f  < < . . . < \+  < 0 . We  now  wish  to  show  < X,  < . . . < \ < 0 , 

1— 2— — n—  12  n 

but  first  we  make  some  preliminary  observations.  We  compute 


x.  -x 


-4,VXj 


fi^fx^x.)  = -[2sinh~2(  1 -■--)  + 3sinh  ‘T(-~)]  (x^)  . 
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and  so  by  the  energy  relation  (6.  1),  4>(x,-x^)  is  bounded.  We  now  claim 
lim  <f(x.-x  ) = 0 . To  see  that,  we  give  without  proof  the  following  easy 

t -.or  ^ 

estimate,  true  for  unbounded  intervals  I, 


(6.  10) 


4 | f ’ | • | f | . > | f j 2 , for  f < L1  fl  L , 

oo  1 — or  7 


where 


f I = J |f  |dt,  | f 1^  = sup  | f | . 
I t€  I 


We  shall  apply  (6.10)  by  picking  our  interval  I = (t,<*),  f(t)  = ^fx.-x^),  and  then 
letting  t -.  «•  we  see  1 4>(x. -xj  | ^ — 0,  as  was  to  be  shown.  Since 

l 1 oc 

Xi"XJ 

f(t)  -*  0,  t — »,  we  must  have  lim  coth  ( — - — ) = ± 1,  depending  on  whether 

t-.sc 

i > j or  i < j . We  are  now  ready  to  show  X.  < X_  < . . . < X < 0 . 

12  n 

Tirst  we  show  X < 0 . By  (6.6),  (6.2),  and  the  asymptotic  behavior 
x -x 

• 1 j X 

of  y .,  coth  ( — - — -),  e 1,  discussed  above,  it  follows  that 


(6.11)  lim  Z±(x(t),  y(t) ) = T*,  lim  W(x(t),  y(t))  = j Tf  T- 

t—  00  t-.  SC 

exists,  where  T*  equals  respectively  a lower,  upper  triangular  matrix, 
with  the  X^'s  in  the  diagonal,  and  ± i respectively  in  all  the  lower,  upper 
entries  respectively.  Hence 

1 n ° n 

det  lim  W(t)  = ^detT+)(detT')  = j(TT  \+)  • (TT  X+)  = - TT  (X.+)  , 
t-oo  i=l  i=l  1 Z i = l 1 

but  since  W(t)  undergoes  an  isospectral  deformation  in  t,  det  W(0)  = 

1 + - X 

det  W(t)  = det  W(°c)  . On  the  other  hand  W = j{Z  )(Z  ' f ae  , but  since 
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-f  ;JC 

(Z  ) = Z , and  as  a > 0,  W is  a positive  definite  Hermitian  matrix, 

which  implies  det  W(^)  = det  W(0)  > 0 . Thus  by  our  evaluation  of 
det  and  since  X^  < 0,  we  have  shown  X^  < 0 for  all  i . 

We  now  show  X*  < X*+1  for  all  i,  using  an  argument  in  [3].  From 
(6.  9)  we  conclude 

e*n  - e'Xl  + (x  -x,)"  > 2 4>(x  -x  ) , 

n 1 — n 1 

x x x1  - 3 

and  since  0<e  n - e I < Cj  e n = 0(t  ) , the  function 

dj  = x - x,  + c,t  1 > x - x,  > 0,  with  some  constant  c_  > 0,  satisfies 
n 1 2 n 1 2 

dt  > 24>(x  -x,)  > 24K  > 0 for  t > t , 

— n 1 — — 0 

where  we  have  made  use  of  4>'(x)  < 0,  x<j>(x)  > 0 . Thus  4j  is  monotonic 

• + + 

increasing  and  ^ (°°)  = X - Xj  > 0 . If  44°°)  = 0,  then  ^(t)  < 0,  and 
since  4>  > 0,  4*  would  be  bounded.  This  implies  4>(^)  is  bounded  away 

from  zero,  which  in  turn  implies  4j(t)  is  unbounded.  This  contradiction 

implies  0 < 4j(°°)  = X^  - X*,  and  since  X*  < X + < . . . < Xn  , there  must 

■f  + + + 

exist  an  s such  that  X < X , . From  this  it  will  follow  that  X,  < X , 

s s+1  1 s 

X+  < X+,  and  proceeding  inductively  we  can  conclude  that  all  the  \.  ' s 
s n 1 

are  different.  We  show  X^  < X^  assuming  s > 1,  the  other  case  being  proven 

in  the  same  way.  Since  X^+  > X* , j > s,  x^  - Xs  > c ^t,  c^  > 0,  t > tQ,  j > s , 

we  have  by  (6.  9)  and  X*  < 0,  that 

.2  _3 

— — (x  -x  ) = ^ 4>(x  -x  )+v  4>(x  -x  ) - 0(t  ) 

dr  j<s  3 J j>i  } 

> 2<t>(xs-x1)  - 0(t‘3)  . 
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Thus  we  are  in  precisely  the  same  position  with  x -x,  as  we  were  in 

s 1 

before  with  x -x  and  so  conclude  X > X.  . We  have  thus  shown 
n 1 si 

+ + + 

X.  < X,  . . . < X <0,  and  similarly  for  t -*■  -°o  i.  e.  , if  x.(-°°)  = X . 
l 2 n ’ l i ’ 

we  find  X.  > x"  > . . . > x"  > 0 . 

1 2 n 


(6.  12) 


Using  arguments  similar  to  those  of  Lemma  2.  1,  we  easily  show 
xi< t)  = X*t  + p*  + G(t  t -*  ± 00,  for  all  i 
y.(t)  = x*  + Oft"2)  , 


and  so  we  must  prove  = -X."  . We  first  prove  this  fact  for  all  |x+| 

i l i 

very  large.  In  that  case  by  (6.  11),  the  spectrum  of  2W  is  'relatively'  close  to 

+ 2 +2  +2 

(Xj)  , (X2)  , ...  ( X ) , and  by  the  implicit  function  theorem  there  exists 
a C'  diffeomorphism  t,  completely  determined  by  (6.11)  and  defined  by 

T((X1V,(X^) (Xn)2>=  ( p2»  • • • » Pn)>  pi  > pz  >»•••>  > pn>  where 

(p^,  p^,  . . . , p^)  equals  the  spectrum  of  W . For  the  domain  N of  t , we 

fix  a point  defined  by  (Xp2  = v.  = (i+l)M,  i = 1,2 M sufficiently 

large  and  positive,  and  let  N = {((X^)2,  . . . , (X^)2  ) | hx^)2  - vj  < p(M)}  , 

P(M)  < 1,  and  being  picked  sufficiently  small  to  ensure  t is  a diffeomorphism. 

On  the  other  hand,  in  the  above  neighborhood,  since  the  spectrum  of  W 

- 2 

is  very  large,  by  ( 6.  11)  in  the  case  of  X , we  know  that  the  (X.  ) , 

i = 1,  . . . , n,  must  be  very  large  and  'relatively1  close  to  the  spectrum  of  W.  In 

fact,  using  the  implicit  function  theorem  once  again,  and  (6.  11),  in  the 

case  of  X , noting  that  the  form  of  (6.  11)  guarantees  that  the  spectrum  of 

W is  unchanged  upon  the  inversion  X.f  — -X.  , i = 1,  2 n.  we  must 

i l ’ ’ ’ ’ ’ 

- 2 - 2 - 2 

have  TflX.  ) , (X.)  , . . . , (X  ) ) - (p.,  . . . , p ) . While  by  the  uniqueness 
i t.  n i n 


-SO- 


4-2  - 2 

clause  in  the  implicit  function  theorem,  this  implies  (X.  ) = (X.  ) , 

and  hence  X,  = -X,  . 

1 l 

f 

We>now  wish  to  prove  X^  = -X.  for  all  (x,  y)  c S2  . For  that  we 
observe  that  (6.  12),  as  in  Lemma  2.  1,  enables  you  to  construct  a canonical 
diffeomorphism  of  S2,  4>  + , defined  by 

4>+(x,  y)  = (X+,  p+),  X+  = (Xj+,  . . . , X*),  p = (Pjf,  • • • , P*)  , 

and  similarly  one  defines  4>(x,  y)  = (X  , p ) . We  can  then  define  b = <b  ° <b  , 
a canonical  diffeomorphism  of  £2,  by  <t(X  ,p  ) = (X  , p + ),  which,  in 
particular,  is  continuous.  Now  by  (6.11), 


(6.13)  W(  ±»r)  r i-  Tf(±  or)  T‘(±oc)  = i 


and  so  the  coefficients  a^  of  the  characteristic  polynomial,  det(zI-W(±sc)) 
n k t- 

= y z a (X*)  depend  on  X or  X,  only.  Thus  since  W(+«>)  ~ Wf-*1)  , 
k=l  k + ' j 

the  map  X -*  X is  given  by  the  n algebraically  independent  relations 
a (X  ) = a.  (X  ),  k = 1,  . . . , n . On  the  other  hand,  we  have  shown  -X  = X^  , 
for  all  k in  some  neighborhood,  and  thus  by  the  continuity  of  4>+,  and  the 

— -f- 

identity  theorem  for  analytic  functions,  -X^  = X^,  k = 1,  . . . , n, 
for  all  S2  . 

To  prove  (6.  8),  we  observe  once  again,  applying  an  argument  of 


Moser' sf  15],  that  the  maps  (x,  y)  — (p*,  x±)  are  canonical,  and  hence 
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II  II 

V dfl+  ~ dx,  = V dx  ^ dy  = y dp  ~ dx  , 

1=1  1 1 1=1  ‘ 1 1=1  1 1 

which  implies,  along  with  the  self  explanatory  definition  of  u , 


du  = d(^  (P*  + pj  )dXi)  = 0,  X.  - X^  - -X.  . 

By  the  Poincare  Lemma  for  the  convex  domain  n,  we  find  u = dS(X)  , 
and  so 

9S(X) 


Pi  + Pi  ^ ax 


i 


We  first  show  S = £ S ( X ),  and  then  we  evaluate  S explicitly. 
1=1 

To  see  S has  the  above  decomposition,  we  note  that  the  one  particle 


system  with 


1 2 x 

H(x,  y)  = - y + e , 


has  for  its  most  general  solution 
(6.14) 


x = log[2X2  sech2(\t+ 6)],  X >0  , 


and  hence 
(6.16) 


x(t)  = ? 2Xt  + p*  + 0(t_1),  t - ± » , 


where  p*  = (log(2X2)  t 26  + 21og2)  . Thus  if  we  have  the  n-particle 
system  with 


(6.  16) 


we  would  have 


v 1 2 x 

«■  E <1  yt  ♦«  *>  > 

i=l  L 
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(6.17) 


x - x = 2 log 

J i 


cosh(\,t  + 6^) 
\,  cosh  ( \,t  + 6 .) 


where  the  notation  obviously  comes  from  (6.  14).  Now  fixing  > x2 

>0  if  6 - 6 , = m > 0,  sufficiently  large,  for  all  i,  we 

n ’ i-1  i 

conclude  from  (6.17)  that  x.  ^(t)  > x^(t)  f°r  aH  N 1 • If  we  let  m "*  r > 

or  what  comes  to  the  same  thing  p|+1  - P^  00  f°r  aH  I then 

inf  (x  -x  ) -*  oo  . We  are  now  in  a position  to  perform  a scattering 
i,  t i+1  1 

experiment  with  the  system  H . Namely,  to  evaluate  S,  since  it  only 
depends  on  the  terminal  velocities,  we  may  choose  a set  of  scattering 
data  for  a fixed  set  of  final  velocities,  i.  e.  , 

x (t)  = -2y  + p*  + 0(t_1),  t - * , 

where  B+  - pf  = m.  m large.  Then  by  the  above  considerations,  the 
Ki+1  i ’ 

system  H decouples  as  m-»oo,  hence 
a 


(6.18) 


= P*  + P“  = nogfZ^)2  + 0(m"1)  . 


But  since  -r^~  (\  , . . . , X ) does  not  depend  on  m,  we  must  have  p.  + p.  : 
d \ 1 n 

2 1 

2 1og(2\.)  . Lifting  the  requirement  that  a = 1,  and  returning  to  the  nota- 

+ - 2 

tion  of  the  theorem,  we  find  p.  + p{  = 2log(o\.)  . 

Corollary  11.  2.  The  explicit  solution  of  the  equations  of  motion  for  a t 0 
are  given  by 


X -1  -1 

(6.  19)  e ~ (la  pp  )’  , 


p = Aj  eAt  A ^ e At  , 


-S3- 


where  ~ indicates  similarity  equivalence.  Here  the  A's  are  constant 

i 

matrices  described  in  the  proof,  and  in  case  a > 0,  we  have 

(6. 2 0)  2 A = diag(-X  -X  ...,  -X  ) 

1 2’  n ’ 

with  X.  being  the  same  as  in  (6.7).  Thus  the  symmetric  polynomials  of 

x ± I X't 

e i,  for  a > 0,  are  rational  in  e 2 J,  j = 1,  . . . , n . Moreover,  if 

we  define  the  2n  X2n  matrices  J,  K by 


— 

0 , I 

K 

y 

0 

iw-  ic_ 

A 

, K = 

0 

y 

K 

where  C =1(1-6^),  and  W(x,  y),  K(x)  as  defined  in  (6.  2),  (6.  3),  then 

(6.22)  Dh  1 = [K,J]  , 

a 

i.e.,  J is  isospectral.  For  a > 0,  the  spectrum  of  J is  ± j X.  , 

i = 1,  . . . , n,  but  in  any  case,  the  spectrum  occurs  in  ± pairs.  Finally  A 

satisfies  the  quadratic  matrix  equations, 

(6.  23)  W ± CA  = 2A2  , 

where 

W = V±*  WV±,  C = V±*CV±,  respectively  , 
for  some  nonsingular  n by  n matrices  V . 

X 

Proof.  We  first  compute  the  time  dependence  of  e , i.e.  (6.19),  and  for 

that  we  go  over  to  (q.,  . . . , q , p.,  . . . , p ) coordinates,  as  in  the  proof 

I n i n ’ 

of  Theorem  11,  using  the  same  notation  as  in  that  proof.  We  define  the 
matrices  q,  r*  by 
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(6.  24) 


q = U_1QU,  r±  = U_1R*  U , 


where  U is  a matrix  satisfying  U(0)  = 1,  U = AU,  A as  in  (6.  3).  This 

matrix  q is  not  to  be  confused  with  the  coordinates  (q, , . . . q ) . The 

1 n 

equations  (6.  4)  are  respectively  equivalent  to  (see  remark  1 of  Section  2), 


(6.  25) 


• - + + 
q = qr  = r q,  since  qr  = r q , 


r = -aq  . 


Since  R+  - R-  = C,  [C,  A]  = 0,  U = AU,  it  follows  that  U^CU  = C,  and 


so  we  have 


r - r"  = C . 


That  W is  isospectral  is  equivalent  to  (see  remark  1 of  Section  2) 

I r+  r"  + o-q  = WQ  , 


where  the  subscript  shall  now,  and  for  the  remainder  of  this  section,  indicate 
evaluation  at  t = 0 . By  the  above  we  have 

1 


(6. 26) 


-(r  + C)r  - r = WQ,  i.  e. 


i,r-,2-r  = W0-tcr' 


Defining  p as  the  solution  of  the  linear  differential  equation, 


(6.27)  r = -2 p p *,  p(0)  = I,  (hence  p(0)  = - j r ' = Z^)  , 


and  substituting  this  into  (6.26),  we  compute 


j(-2pp  Sf-Zpp'1)  . ( -2  p p" 1 + 2pp'Ipp'1)  = WQ  - — C(  -2  pp1)  , 


which  implies 


5 =tw0p,tcp 


Defining  the  n by  2n  matrix  P=(£),  we  conclude 


(6.  28) 


V 


P = J P , hence  P = e P„  , 

U 0 7 


with  J defined  in  the  statement  of  this  corollary. 

Recalling  Q = diagfq^  • . . , = diag(exl,  . . . , exn)  = eX,  see  (6.  3), 

we  compute,  using  (6.24),  (6.25),  and  (6.27), 

(6.29)  U e U = U QU  = q = -a  ^r  = -a  \-2pp  V = 2a  *(pp  ^)’,  hence 

U'1eX(J  = 2«'1(pp'V  . 


Assuming  for  the  moment  the  statement  of  the  corollary  concerning  the  spectrum 
of  J,  we  may  write,  using  (6.28), 


E e—  1 


(E 


po> 


9 


with  A = diag(A,  -A),  A an  n Xn  matrix  which  takes,  for  a > 0,  the  form 

given  in  (6.  20),  and  E a matrix  which  block  diagonalizes  JQ  . This  expression 

for  P,  in  conjunction  with  the  above  expression  for  UeXU*1,  yields  (6.  19). 

We  now  analyze  the  spectrum  of  J,  but  first  we  must  verify  (6.  22), 

i.  e.  J = [K,  J]  . We  compute 

a 
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0 


0 


dh  j = 

a 


2 °H 


while  by  (6.21), 


[K,W],  [K,q  [K,W], 


and  thus  by  Theorem  11,  (6.22)  is  verified. 

We  now  assume  a > 0,  returning  to  (x,  y)  coordinates,  and  proceed 
to  prove  that  the  spectrum  of  I (which  by  (6.  22)  is  constant  along  an 
orbit)  is  ± j i = 1,  . . . , n . Now,  since  J(t)  is  an  isospectral  matrix, 

it  is  sufficient  to  prove 


J(  op)  = lim  J(t)  = 

t-»  00 


5-tY, 


(see  (6.13)),  has  the  above  spectrum.  Defining  the  matrices  A,  X,  by 


Ajk  = 


i if  k < j 


0 if  k > ) 


, \ = diagfiy  V,,  . . . , M , 


(6.  13)  implies 


T+T_  = ( \f A)( X+A  ),  C = A - A 
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Thus  to  prove  that  the  spectrum  of  J(oo)  js  + L \ j _ i n . 

' *2  i.9  1 ~~  A 9 " " * * *■ 

sufficient  to  prove  the  seemingly  stronger,  purely  algebraic  statement; 
Lemma  11.  3.  If 


(6.  30)  M (X)  = 
e 


eA)(\+eA  ), 


I 

1 * » 
- e(A  - a ) 


where  X = diag( Xj,  . . . , X^,  the  spectrum  of  M£(X)  is  independent  of  e, 
and  in  particular  the  spectrum  of  MQ(X)  agrees  with  the  spectrum  of 
Mj(X)  = J(»)  . 


ProoL  r°r  the  Proof,  which  need  only  be  given  for  the  case  X < X < 

1 2 

' * ‘ < Xn  < 0,  we  emPloy  the  asymptotic  description  of  the  orbits  of  H 

< 

formula  (6.7).  First  we  make  some  preliminary  observations.  Defining 


X = X(e) 


0 

el 


we  compute 


(6.  31) 


e BM  ( X)B  = M (X'),  i.  e. , eM.(X)  ~ M (X')  . 

e 1 p 


For  the  rest  of  the  proof  of  Lemma  11.  3,  we  fix  X',  Xj  < X'  < . . . 

< xn  < 0 . We  then  pick  6j  = S^X1)  so  that 

(6.32)  |Me(X‘)  - M0(X’)I  < J&2,  if  UI<6J<62<I  , 

with  &2  = 62(X ')  to  be  determined  in  the  course  of  the  discussion,  and 
n 

U|  = sup  (^  I A | ) . We  now  fix  e for  the  remainder  of  this  lemma. 


1=1, . . . , n j = l 


e only  having  the  property  |e  I < 6 
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Observe  that  for  any  matrix  of  the  form  M^X')  + S,  S a small 


matrix  with  |S|  < 6 , we  may  find  a matrix  D such  that 


(6.  33) 


D-  (M0,v,  fS)D  = dl.g(Ml+ 


o 


± l 


> 9-  i = J + °(62))>  i = 1,  • • • , n 


In  addition,  D may  be  taken  to  be  of  the  form 


r 


D 


rd,,. 

d,  i 

i , i 

a* 

12 

d 

22 

Lx-  L v 
'2  K ’ 2 v 

+ 0(62)  , 

(6.  34)  / and  so 


D-‘  = 


I1* 

-(X1) 

I1* 

(V) 

-1 


-1 


+ 


Note  that  all  the  O's  which  appear  now,  and  for  the  rest  of  the  discussion, 
depend  only  on  X*  . 

We  proceed  with  the  proof  of  Lemma  11.  3,  having  fixed  X‘,  and  then 
e,  and  therefore  X = e X',  as  previously  discussed.  We  first  consider  a 
particular  orbit  with  asymptotic  behavior  given  by  (6.7),  and  any  choice  of 
P = (P|.  • • • P ) • We  then  pick  the  origin  of  the  orbit,  i.  e.  , replacing  t by 
t + t , with  t chosen  so  large  that 


(6.  35) 


6 

"J0(x(t),  y(t))  = Mj(X)  + Sj,  |Sj|  <y 
|^Z‘(x(t),  y(t))  = (X  +A  :)  + S2>  |S2I  < 1 


-59- 


Since  by  (6.  7), 


lim  J(x(t),  y(t))  = M (\),  lim  Z"(x(t),  y(t))  = \ +A  , 

t -*  or  t -*■  00 


such  a t can  always  be  chosen. 

With  the  above  chosen  orbit,  we  rescale  time,  t — e ^t  = t ' , and 
so  in  (6.7),  the  terminal  velocities  X are  replaced  by  eX  = X'  . In  what 
follows,  we  shall  indicate  the  transformed  quantities  by  primed  letters. 


We  compute 


(6.  36) 


( P.  ) = P - BP  = P1  = (P.),  B = 
P P 


[o,  el]  ’ 


p=  I0P-  eBP  = (eBJ0B  )(BP)  , 


which  may  be  conveniently  written  as 


(6.  37) 


3F  = Jo  P’’  Jo  = eBIoB'1 


Using  (6.  35),  (6.  31),  we  compute 


VQ  = eBIo8"1  = eB(Mj(X)  + S^B'1 


= M (X1)  + eBS  B = M ( X')  + S,  + S,  , 
e 1 0 3 4 


S3  = Me(X')  - Mq(X'),  S4  = eBSjB 


Since  by  (6.  32),  I e I < 1,  I S4  I < I B I • I S{  | • | eB-1 1 < | S I < j 6 , , 


and  so,  again  by  (6.  32),  we  find 


(6.  38)  l'Q  = M Q(  V)  + S,  |s|  < . 

We  are  now  in  a position  to  use  (6.  33),  as  (6.  37)  implies,  with  D de- 
fined in  (6.  33), 

J' t , j 

P'  = e 0 Pq  = D {exp[(D  J^D)t]>  (D~  P^)  , 


and  so  by  (6.  38)  applied  to  (6.  33),  (6.  34),  we  find 


We  shall  new  come  to  the  determination  of  6^  = 6^(\'),  and  then  we 
will  finish  the  proof  of  this  lemma.  By  (6.  36),  (6.  27), 


I 


I --(X.A 


+ S. 


Since  by  (6.  32),  (6.  3S),  | e | < 6^,  | | < 1,  we  conclude  from  the 

above 

(6.41)  P,  = I + S , with  I S | < <r  ( 6 ) , 

ID  D L 

*(&-,)  = jlv  I"1  6z(nfl)  f 0(62)(1  |XJ|  4 j 6z(n  +1))  . 

We  now  choose  &2  = 6 ^ ( V ) , so  that  0 < &2  < 1,  o-(62)  < 1,  and  in  addition 
&2  is  in  the  domain  of  validity  of  the  estimates  in  (6.  33),  (6.  34)  . The 
remark  concerning  the  V dependence  of  O ensures  &2  = &2(\')  • 

-1  X 

We  now  compute  the  positive  spectrum  of  . Recalling  q = U e U , 
(see  (6.29),  and  substituting  t — t* , \ in  the  asymptotic  description 

of  the  orbit,  (6.7),  we  conclude,  using  (6.  24) 

U(t)  = U(°°)  + OUt')"1),  lim  U(t')  = U(°o)  , 

t 1 -*•  00 


q.  = exi  = e0i  e V(1  + 0((t')-1)) 


qt  = eXi  = e0i  e ^(l  f 0((t1)'1) 
i = 1, . . . , n,  and  so 

(6.42)  q^q  = [U_1(*)  (diag(\j,  . . . , V ))  U(=o)]  • (1  + 0((t')_1))  , 


where  these  O's  depend  on  our  fixed  initial  data.  On  the  other  hand,  using 
(6.40),  and  the  nonsingularity  of  the  matrices  d^,  p implicit  in  (6.  34),  (6.  41) 
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respectively,  and  since  by  (6.25),  (6.27),  q q = r = -2pp  , we  conclude 

-1.  . -1  , + -1 

lim  q q = lim  (-2pp  ) = -2d ^ p d . 

t'-»00  t'-^00 

This  expression  coupled  with  (6.42)  yields 

U 1(<»)(diag(\|,  . . . , \^)>  U(»)  = -2  du  p f d^1  . 

We  thus  find  that  the  positive  eigenvalues  of  J^,  p * = - j Xj  , 
i = 1,  . . . , n,  with  no  error  term  0(6^),  as  in  (6.  33).  Thus  p , , 

i = 1,  . . . , n,  is  contained  in  the  spectrum  of,  (see  6.  37), 

J*  = eBJ  B'1  ~ eBJMB'1,  JM  = lim  J(x(t),  y(t))  , 

U U t-*oo 

where  ~ denotes  similarity  equivalence,  since  J(x(t),  y(t))  undergoes 
an  isospectral  deformation  by  (6.  22).  By  (6.7),  K00)  = M^(X),  and  thus  by 

(6.  31), 

JL  ~ eBMdMB'1  * M (V)  . 

0 1 £ 

We  have  thus  shown  for  arbitrary  e satisfying  |e|  < 6^  = 6 ^ ( \* ) , that  the 

positive  spectrum  of  M (V)  contains  - — X'  i = 1,  . . . , n . By  the  same 

e c * 

argument  for  t’  -*  -<»,  after  having  first  shifted  our  origin  appropriately, 

we  also  conclude  that  -j  \'f  i = 1,  . . . , n is  contained  in  the  spectrum  of 

M (V)  . Hence,  M (V)  has  full  spectrum  ± 7 X'  i = 1,  • • • , n,  at  least 

e e ci 

for  |el  < 6^  = 6j(V)  . We  therefore  have  proven  Lemma  11.  3 for  e small,  but 

since  M (X)  is  a polynomial  in  e,  the  proof  is  finished, 
e 

To  complete  the  proof  of  Lemma  11.  2,  we  must  verify  (6.  2 3).  To 
accomplish  this  end,  we  shall  construct  an  analytic  invertible  matrix 
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function  of  (x,  y),  which  diagonalizes  J,  namely 


(6.43)  E = E(x,  y)  = 

possessing  the  property  that  e..,  i,  j = 1,2,  is  invertible.  Assuming  for 
the  moment  the  construction  of  such  an  E,  by  (6.  21),  (6.  20),  and  the 
spectral  resolution  of  J,  we  have 


which  implies 


which  clearly  implies  (6.23),  with  V,  = e^,  V = . Note  by  (6.45) 

and  the  nonsingularity  of  A,  that  once  we  have  constructed  E(x,  y) 
satisfying  (6.44),  the  invertibility  of  e^,  e implies  the  invertibility  of 

eu,  M = 1,2  . 

First  we  shall  construct  E(x,  y)  satisfying  (6.44),  and  then  proceed 
to  show  the  invertibility  of  its  blocks.  The  construction  shall  proceed  in  two 
steps,  i.e.,  we  shall  take  E of  the  form  E^(x,  y)  E,(\)  . To  construct 
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E we  first  observe  that  using  the  arguments  of  remark  2 of  Section  2,  we 
may  algebraically  construct  a global  unitary  matrix  U (x,  y)  such  that 


(6.46)  U \x,  y)  W(x,y)U(x,y)  = j(X+A)(X+A  ),  C U - U C - C, 

— 1 * 

C = C + il,  with  X.,  A defined  as  in  Lemma  11.  3,  y(X+A)(X+A  ) being  the 

limiting  value  of  W(x,  y)  along  an  orbit  with  initial  data  (x,  y).  It  follows 

by  (6.21),  (6.46),  that  if  we  define  the  2n  X 2n  unitary  matrix 


El  = 


U,  0 

L°  » u 


,-i 


, then  Ej  JEj  = Mj(\)  , 


Mj(M  defined  in  (6.  30)  . Moreover,  we  may  uniquely  construct  an  invertible 


analytic  matrix  = E^(\),  so  that 


e2  Mi(X)E2  = 


A , 0 

0 , -A 


, A as  in  (6. 20)  , 


with  E^  normalized  in  the  following  way*.  Pick  a X'  and  fix  it,  defining 
X = e * X',  but  unlike  the  situation  in  Lemma  11.  3,  we  shall  think  of  e — 0 . 
By  (6.  31),  (6.  33),  we  have 

eBM,(X)B_l  = M (X1)  = M IX')  + O(e) 

1 e u 

= D A'  D , with  A'  = 

with  O depending  on  X',  and  where  we  have  employed  Lemma  11.  3 to  con- 
clude A ' = - — X'  . 

From  this  we  infer 
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(6.  17) 


(B'1  D)"1  Mj(X)  (B_1D)  = A , A=yX, 


where  by  (6.  34),  and  B - 


I , 0 

0,  el 


, we  conclude 


r 

-1 

~I,  0 “ 

f 

B D = 

0,  e_1I_ 

\ 

I , I 

-lr  X'.  IrX' 


f O(e) 


"i,i" 

~0(e),  O(e) 

i „ i , 

+ 

-1  -1 

- zr  ~ X 

L 2 ’ 2 J 

e O(e),  e 0(e) 

O depending  on  V,  and  similarly  from  (6.  34)  we  find 


-1  -1 


1 

, -i 

1 „ 

21’ 

"O(e), 

eO(e) 

B = 

|> 

X-1 

w 

+ 

_0(e), 

eO(  e) 

The  above  estimates  allow  us  to  normalise  E^,  so  that,  again  keeping  V 


fixed,  and  thus  fixing  the  behavior  of  O(e)  as  e -*•  0 , 


0 as  O(e)  , 


or  equivalently 


0 

-1 
e I 


0 as  e -»  0 as  O(e)  • 


We  thus  have  constructed  the  invertible  analytic  matric  function 

E(x,  y)  = E (x,  y)  • E7(X),  and  so  we  have  from  (6.  28),  the  time  dependence 

of  the  H flow  (6. 1)  given  by 
a 
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( P)  = P = E ■ {exp[(E_1  JQ E)t ] } ■ (E*1Pq)  , 


with  E depending  on  the  initial  data)  , 

I 1 


Ee-‘<E'V  V 


\ zo 


eir  ei2 

’eAt,  0 ’ 

"Pl" 

X 

e21’  e22 

0 , e'At 

_P2_ 

y 

A 

, by  (6.44),  (6.27)  , 


,-l 


= E PQ,  i.e., 


(6.48) 


At  -At 

P=eil6  Pl+ei26  P2  • 


Note  that  the  invertibility  and  analyticity  of  E implies  that  e.^  , 
i,  j = 1,2,  Pj,  p^  are  analytic  functions  of  the  initial  data.  By  arguments 
similar  to  the  ones  given  above,  and  in  Lemma  11.  3,  we  can  easily  show 
that  e^,  Pj  are  invertible  on  some  open  set  in  (x,  y)  space,  and  similarly 
for  e , p^  . We  thus  conclude  the  proof  of  (6.23),  and  therefore  of 
Lemma  11.  2,  by  proving  the  following  statement; 

Corollary  11.  4.  Under  the  hypothesis  of  Theorem  11,  we  have 

n n 

+ 


(6.49) 


(a)  Y 3.  = V x.(0)  - 2 log(det(e  p )) 

i . l 11  i 

1=1  1=1 

n n 

(b)  V P~  = }]  xi(0)  ' 1 lo<3(det(elzP2)) 


a > 0 , 


with  6^,  e^,  Pj,  p^  as  in  (6.48),  p , p"  as  in  (6.7). 

Remark.  This  easily  implies  e * exists  for  i,  j = 1,2,  by  the  aforementioned 
analyticity  of  e , p , p , and  of  p , in  initial  data. 

1J  1 £ 1 
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Proof  of  Corollary  11.  4.  We  shall  just  prove  (a),  (b)  follows  similarly. 


Clearly  it  suffices  to  establish  (6.49a)  in  some  open  set  where  e^p^ 

. - 1 r 

is  nonsingular,  and  then  analytically  continue.  Since  pp  = , 


p(0)  - I,  we  compute,  using  (6.  3), 

d log  (detp)  1 t 1 t 1 v 

-± = - j tr(r  ) = - j trR  *-J  } Y,  ■ 

i=l 

Hence  log  detp  4-  j V x.  = constant,  and  thus 

\ I'  xi(0)  = lim  (log(detp)  + l-  Yj  x .)  = 
t -*  on 

n - Y 

lim  {logfdetfe^))  + log(rr"  ! e— — +■ 

t C30 

1 n 1 + - 1 ^ 

(2  L Y +Z  l Pi  + °(t  )}  ’ 


1 n 

7 V *■ 

z iti  1 


where  we  have  used  (6.48),  (6.7),  from  which  we  conclude 

i-  Y Y°>  = log(det(euPl))  + j Y?{  > 

and  the  corollary  is  proven. 


-68- 


REFERENCES 


F.  Calogero  and  C.  Marchioro:  Exact  solution  of  a one-dimensional 
three-body  scattering  problem  with  two-body  and/or  three-body  in- 
verse-square potential,  J.  Math.  Phys.  15(1974),  1425-1430. 

B.  Sutherland:  Exact  results  for  a quantum  many-body  problem  in  one 
dimension,  II  Phys.  Rev.  A 5(1972),  1372-1376. 

J.  Moser:  Three  Integrable  Hamiltonian  Systems  Connected  with 
Isospectral  Deformations,  Adv.  in  Math.  vol.  16,  No.  2,  May  (1975). 

P.  D.  Lax;  Integrals  of  nonlinear  equations  of  evolution  and  solitary 
waves,  Comm.  Pure  Appl.  Math.  21,  (1968). 

H.  Flaschka;  The  Toda  Lattice  I,  Phys.  Rev.  B (1974)  , 1924-1925. 

Arnold  and  Avez,  A.  : Problemes  erogodiques  de  la  mecanique  classiaue. 
Paris  1967.  English  translation;  Ergodic  problems  of  classical 
mechanics.  New  York  1968.  See  Appendix. 

G.  Gallavotti  and  Marchioro:  On  the  Calculation  of  an  Integral, 

J.  Math.  Anal.  Appl.  44(1973),  661-675. 

A.  Weinstein;  Normal  Modes  for  Nonlinear  Hamiltonian  Systems, 
Inventiones  Math.  20,  47-57,  (197  3). 

K.  Sawada  and  T.  Kotera;  Integrability  and  a Solution  for  One 
Dimensional  N- Particle  Systems  with  Inversely  Quadratic  Pair 
Potentials,  preprint,  Tokyo  University  of  Ed.  TUETP-75-10,  See 
Appendix 


-69- 


[101  F.  Calogero:  Solution  of  the  One -Dimensional  N-Body  problem 

with  Quadratic  and/or  inversely  quadratic  pair  potentials,  Journal 
of  Math.  Phys.  12,  (1971),  419-436,  esp.  p.  427. 

[11]  M.  A.  Olshanetzky  and  A.  M.  Perelomov:  Explicit  Solution  of  the 
Calogero  Model  in  the  Classical  Case  and  Geodesic  Flows  on 
Symmetric  Spaces  of  Zero  Curvature,  Lettere  A1  Nuovo  Cimento  16, 

333,  1976. 

[12]  M.  Adler;  A New  Integrable  System  and  a Conjecture  by  Calogero, 
Nov.  (197  5)  unpublished. 

[13]  M.  Adler,  Some  Finite  Dimensional  Integrable  Systems,  preprint 
6176,  to  appear  in  proceeding  of  conference  on  Theory  and  Application 
of  Solitons,  held  1/76. 

[14]  M.  Adler,  Some  Finite  Dimensional  Integrable  Systems,  doctoral 
dissertation,  Courant  Institute,  10/76. 

[16]  J.  Moser;  The  Scattering  Problem  For  Some  Particle  Systems  On  The 
Line,  to  appear  in  the  proceedings  of  the  Conference  on  Dynamical 
Systems,  held  in  Rio  de  Janeiro,  summer  of  1976. 


-7  0- 


